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FUNÇÕES DELTA NO UNIVERSO 
DA ANALISE NAO-STANDARD 


J. Sousa Pinto * 


1. INTRODUÇÃO 


Dois dos mais bem conhecidos e significativos exemplos de distribuições de 
Schwartz são, respectivamente, a “função” salto unitário de Heaviside, H, e 
a “função” delta, ô, que “funciona como derivada” de H num sentido que se 
procurará explicitar. Embora Heaviside (1850-1925) tenha usado ele próprio 
ambas as funções H e ô nos seus trabalhos, ' a verdade é que as proprieda- 
des mais características atribuídas à função delta bem como a notação ó(x), 
hoje universalmente adoptada, só apareceram explicitamente nos anos 20 
deste século num texto de Mecânica Quântica de Paul M. Dirac [I]; por esta 
razão, normalmente, chama-se-lhe função delta de Dirac. O nascimento da 
função delta, no entanto, parece poder situar-se anteriormente tanto a Dirac 
como a Heaviside. No seu livro sobre a pré-história da Teoria das Distribui- 
ções, Lútzen [2] sugere o aparecimento da função delta em 1822, num texto 
de Fourier (1768-1830) intitulado “Théorie Analitique de la Chaleur”. 

Em livros de texto elementares (de física e engenharia) é frequente apre- 
sentar-se a função delta “definida pontualmente” da seguinte forma 


sl? se 140 
WE, se x=0, (1) 
assumindo-se adicionalmente a condição 
fostja=1 (2) 


Nenhuma função, no sentido rigoroso do termo que lhe é dado em Análise 
Matemática, pode satisfazer este comportamento. No entanto, embora fa- 
zendo uso de manipulações puramente formais, (1) e (2) podem justificar-se 
da seguinte forma: seja 4 : IR > IR a função definida por 


(o) se «<0 
H(x)=te se x=0 = 
1 se x>0, (3) 


onde «e(0.1). A função derivada de H. (x) é função nula para todo o 140 
e não é definida para x=0. No entanto, a derivada de H (x) na origem 
existe, sendo 


HA)=H(0) 


H:(0)= lim 
a x 


=+eo, 


* Departamento de Matemática da Universidade de Aveiro 

! Os símbolos H e 5 não designam funções no sentido rigoroso da palavra. No entanto, 
porque são frequentemente referidos como tal na literatura técnica, adoptaremos aqui 
a mesma terminologia (embora isso constitua, claramente, um abuso de linguagem). 


Sendo 5(x), um objecto matemático ao qual será atribuido o papel operacio- 
nal de “derivada” de 4 (x) então os resultados obtidos acima servem de 
motivação para a anómala definição (1). Com a mesma lógica (ou falta 
dela!) com que se considera ô a “derivada” de H, pode entender-se que 
H, é uma “primitiva” de ô; consequentemente, uma integração formal da 
função delta, originaria o resultado seguinte 


| atdr= tipo AH (lj = 


que serve de justificação para (2). 

Embora formais e sem sentido matemático, a verdade é que as funções H e 
à gozando, de algum modo, as propriedades acima referidas, têm enorme 
importância nas aplicações. Assim vale a pena abordar o seu estudo e pro- 
curar o que se esconde por detrás das suas (aparentemente) contraditórias 
propriedades. Procuraremos fazê-lo aqui sem sair de um contexto relativa- 
mente elementar de Análise Matemática acessível ao leitor não 
especializado ?. Como objecto matemático não pode aceitar-se a definição 
dada acima e, por este facto, no contexto da Análise Matemática (clássica) 
só de uma forma meramente simbólica se pode empregar a palavra função 
quando esta se refere ao símbolo ô. No entanto, se em vez do universo da 
Análise Standard nos colocarmos no universo da Análise Não-Standard, à 
pode ser representada (embora de forma não única) por uma função no 
sentido rigoroso do termo, podendo mesmo requerer-se que tal função seja 
indefinidamente diferenciável. 


2. A FUNÇÃO DELTA DE DIRAC 


2.1. Significado Operacional 


Recorde-se, antes de mais, que se chama suporte de uma função 
p:lR > ao fecho do conjunto dos pontos onde q não se anula. Sendo p 
um número inteiro não negativo denota-se por C/=C/(IR) o espaço vecto- 
rial de todas as funções continuamente diferenciáveis pelo menos até à or- 
dem p cujo suporte é compacto * O espaço das funções pertencentes a c/ 
qualquer que seja pe N denota-se geralmente por Cj.* 


Considere-se a função /, definida como se indica em (3) onde c designa 
agora um número real qualquer. Sendo «p uma função em C,, o integral 


fo moolidr= [ola (4) 


é bem definido qualquer que seja y. O segundo membro de (4) mostra 
mesmo que o valor daquele integral não depende da constante «e Ik consi- 


Para uma abordagem mais profunda o leitor interessado pode consultar, por exem- 
plo, as referências bibliográficas e, em particular, o livro do Professor Campos Fer- 
reira recentemente editado pela Fundação Calouste Gulbenkian [3]. 

Um compacto de /R é um conjunto limitado e fechado. 

Este espaço tem elementos diferentes da função nula, como se pode provar, por 
exemplo, com a função de Schwartz q:IR- IR definida por elx) =exp( 
hi<te q(x)=0 se |x>1. 


derada em H, Pode então identificar-se toda a família de funções local- 
mente integráveis [H.] ., com a funcional H:Cy > definida por 


<H,9>=Hlo]= [o(vds 


Esta funcional é independente do valor de c razão pela qual se pode deno- 
tar simplesmente por H *. Como já foi referido, por abuso de linguagem, 
fala-se muitas vezes na função H, mas em bom rigor esta designação de 
função só pode ser legitimada se entendida num sentido generalizado. Qual- 
quer dos elementos !, pode ser usado para representar H símbolo este 
que, quando muito, especifica uma classe de equivalência de funções. 

A identificação de uma família de funções com uma funcional é frequente 
em Análise Matemática. Qualquer função localmente integrável f:k>q dá 
origem a uma função generalizada F entendida naquele sentido, isto é, a 
uma funcional F:C, > definida por 


<Fp>=F[o]=[" o(a)f(a)a (5) 


O símbolo F denota também uma classe de equivalência de funções visto 
que se pode substituir f por outra qualquer função f,:IR=> desde que se 
tenha f(x)=f(x) em quase toda a parte (isto é, em todo o IR excepto, 
quando muito, num conjunto de medida nula). Se f for uma função contínua 
pode, sem qualquer perigo de confusão, identificar-se a funcional F com fjá 
que, neste caso, a classe de equivalência F contém apenas um elemento 
com aquela propriedade. 

Suponha-se agora que f pertence a C' e seja y uma função também conti- 
nuamente diferenciável de suporte compacto contido, por exemplo, no inter- 
valo [a.bJc IR; fazendo uma integração por partes e atendendo a que 
p(a)=q(b)=0 vem 


Poridogar= |" rodo aja. (6) 
A função /”, sendo contínua, gera também uma funcional sobre o espaço 
€; pelo processo atrás indicado que, pelas razões então invocadas, se 
pode identificar com 1”. Assim sendo, é natural designar esta funcional por 
derivada da funcional F gerada por fe, com este sentido, representá-la pelo 
símbolo DF. Da equação (6) resulta que a funcional DF:C; >C se pode, 
neste caso, definir para qualquer pec, por 
<DEg>=[" [(e)o(a)dr (7) 
ou, indiferentemente, por 


<DE,p>=[. Flo (Jar (8) 


Escolhendo a equação (8) para definir DF, pode então estender-se o opera- 
dor de derivação D a toda a funcional F gerada por uma função f que seja 


* Este facto, provavelmente, justifica que em muitos livros de texto elementares a de- 
finição da função de Heaviside seja restrita ao conjunto 1R (0). 


apenas localmente integrável; neste caso, fe; portanto o segundo mem- 
bro de (7), pode não fazer sentido, mas o integral do segundo membro de 
(8) é sempre bem definido. Por exemplo, 


[+vx se x20 


Lo se x<o 


é uma função absolutamente contínua, diferenciável no sentido usual ape- 
nas em JRN([0]. A sua derivada clássica é a função 


J+y2vx se x>0 


l 
dugk = 
2 [o se x<0, 


que não é limitada em qualquer vizinhança da origem. No entanto, qualquer 
que seja a função qyeC, as funções v e 4x, verificam a equação (6). 
Então 4x. é um representante da derivada DF da funcional F gerada pela 
função ordinária Ki - 

Considere-se de novo a equação (6) quando se substitui f por um represen- 
tante qualquer 4 da função generalizada H. A função (x) admite como 
derivada ordinária, qualquer que seja ce !k, a função que é nula em todos 


os pontos 140 e não é definida no ponto x=0. Então do lado esquerdo 
da equação (6), para qualquer pec,, obtém-se 


f Le) Cada = lim |" (sto NH (x )de = 0; 


do lado direito de (6), supondo que j tem suporte contido no intervalo [a.b], 
obtém-se 


Il [Eorto]a (jar -[ [ota =[-9"(9], = (0) 


resultado este que, de novo, se verifica ser independente da constante 
celR. A fórmula (6) de integração por partes deixou de verificar-se. No en- 
tanto a definição de derivada de uma funcional sobre c, dada em (8) per- 
mite estender à funcional H a noção de derivada obtendo-se a funcional de- 
rivada DH:C) = € que é definida por 


<DH.q> 


PE a of-o(a)Jdr = (0) 


para qualquer qec,. Esta funcional, que realiza a chamada operação de 
amostragem + (1), é normalmente denotada pelo símbolo 6, tendo-se 
então 


6=DH, 
podendo dizer-se, com este sentido, que a função delta de Dirac é a deri- 


vada (generalizada) da função de salto unitário de Heaviside . Se se inter- 
pretar a expressão 


faceta [=] DH(a)o( ds | 


num sentido meramente simbólico, atribuindo-lhe, por definição, o valor que 


a funcional derivada, DH, assume no ponto yeC,, pode então restaurar-se 
a validade formal da fórmula de integração por partes e escrever 


E (x) p(x)dr =[. Hlo)f=9'()Jdr 
onde o integral no segundo membro é entendido no sentido usual quando 


se substitue H por uma sua representante /H arbitrária. Daqui resulta a ex- 
pressão simbólica 


fo otsótedr = (0). (9) 


Considerando, em particular, a função «(x)=1, V,., obtém-se a formula (2) 
proposta (heuristicamente) por Dirac. Observe-se, finalmente, que à fórmula 
(9) se pode dar um sentido mais geral podendo mesmo ser escrita para 
qualquer função y que seja apenas contínua numa vizinhança da origem. 


2.2. Derivação Generalizada. 


Visto que C; está contido em C; qualquer que seja pe IN,, então tudo o 
que acima se disse permanece válido quando se consideram apenas fun- 
ções em C;. Sem qualquer perda de generalidade, e por uma questão de 
simplicidade, tomar-se-á C, para domínio das funcionais consideradas a 
menos que, explicitamente, se diga o contrário. Tendo (8) como motivação, 
pode agora introduzir-se a seguinte definição: 


Definição 2.1 Dada uma funcional F definida em C, chama-se derivada 
generalizada de F à funcional DF definida por 


<DF,p>=(DF)[p]=F[-q]=<F.-q'> 
qualquer que seja pec; . 
Seja x,:IR > IR a função definida por 
se x>0 


safa 
“lo se xeo 


Visto tratar-se de uma função contínua em !R pode identificar-se, sem qual- 
quer perigo de confusão, à funcional que ela gera em C;. Denotando esta 
funcional pelo mesmo símbolo, tem-se então 


xp(x)dx 
qualquer que seja qec,. A derivada generalizada de x, é a funcional 
Dx,:C; > definida por 


De foJ=<4.-q'>= [" af=o(a)Jde= [o(ajar 


o que mostra que a funcional Dx. coincide com a função salto unitário de 
Heaviside, H. Com este sentido pode escrever-se 


H=Dr. 


10 


A derivada de Dx, é a funcional D'x,:C; > definida por 
D'clo]=<Dx,-p'>=<x,.9">= [xp"(U)de=9(0) 


o que mostra que D'x. coincide com a função delta de Dirac. Com este 
sentido pode também escrever-se 


ô=DH=D'x., 


o que significa que a função delta de Dirac é igual à segunda derivada ge- 
neralizada da função contínua «,. 


Aplicando agora a definição 2.1 à função delta, vem 
Dôlp]=<Dô.p>=<3.-q'>=-q(0) (10) 


qualquer que seja q eC;. A derivada generalizada Dô da funcional ô re- 
presenta-se frequentemente por 5” embora esta notação tenha o inconveni- 
ente de se confundir com a operação usual de derivação. Esta prática está, 
no entanto, muito vulgarizada e, de facto, não apresenta grandes perigos de 
confusão pois fica claro a partir do contexto qual o sentido que deve ser 
dado à operação. Por indução finita pode definir-se a derivada generalizada 
de ordem me IN, da funcional ô obtendo-se 


D"ô[p]=6"[p]=<3(1)" q" >=(1)" q!” (0) (11) 


qualquer que seja pec; A fórmula (11) pode interpretar-se simbolicamente 
da seguinte forma 


fo sr (otdr=(-" 0" (0) 
que é uma fórmula generalizada, também proposta por Dirac. Pode então 
concluir-se que a função delta de Dirac é indefinidamente derivável e ainda 


que a derivada 6” de qualquer ordem me IN, é a derivada (generalizada) 
de ordem m + 2 da função continua x,. 


3. REPRESENTAÇÕES NÃO STANDARD 


3.1. Limites Generalizados e “integração 


Considere-se, por exemplo, a sucessão de funções d,:IR > IR, n=1,2.... de- 
finidas por 
| nt se 0<r<ljn 
d(t)=/2n-n't se Yn<t<2n 
| 0 se te(0.2/n). 


O limite pontual desta sucessão de funções é a função nula em toda a 
parte, ou seja, lim,..d,(1)=0. V,,- Esta sucessão de funções possui, no 
entanto, algumas propriedades interessantes. 


(1) Visto que o integral de d, estendido a IR é igual a 1, para todo neN, 
então 


tim [º d(dt=1 


enquanto que o integral do limite pontual de d, (quando n=><) é igual a 
zero. Este fenómeno, bem conhecido, resulta do facto de os operadores de in- 
tegração e limite não serem, em geral, comutáveis. De qualquer forma não 
deixa de causar alguma perplexidade (ao nível da intuição) o facto de ser dife- 
rente de zero o limite do integral de uma sucessão de funções que tende (pon- 
tualmente) em toda a parte para a função nula. Mas tal perplexidade, se assim 
se pode dizer, é ainda mais justificável quando se considera um limite do tipo 


tim [o a (op(idr 


onde «p designa uma função contínua numa vizinhança da origem. Visto que, 
e(o)- |" atne(iar= |" [o(0)-o(n]a, (rar 
=[" [o(o)-o(n)rrdr+ |” [o(0)-0(0)[2n=nºr)ar 
e, como para n suficientemente grande, se pode supor que q(0)-q(t) é 


uma função contínua no intervalo ]0.2/n[ então, aplicando o Teorema do 
Valor Médio generalizado, vem 


o(o)-["a(op(rar=+[o(0)-p(E,)] +[9(0)-0(1,)] 


onde 0<é,<l/n e I/n<n,<2/n. Fazendo n tender para infinito Sen, 
tendem para O e, portanto, 


tim Pastora = (0), 


o que permite concluir que, embora a sucessão (d,),.. “tenda a desapare- 
cer” quando n tende para infinito, não o faz sem deixar “marcas” em todas 
as funções q contínuas numa vizinhança da origem, tomadas para funções 
teste. Procuremos agora uma interpretação não-standard para estes factos. 

Seja (/,),.« uma sucessão de funções standard cada uma das quais é defi- 
nida e integrável no conjunto 4, c |. Pode então definir-se o “integral da 


função interna F=[(1,), dei] estendido ao conjunto 4=[(4,),.,]= "It como 
sendo o número hiper-real 


“fi rtoa= ( flejar | RI] 


As propriedades do integral (standard) são transferidas para o “integral de 
funções internas; em particular, facilmente se estabelece uma versão 
não standard para os Teoremas da Média. Se 4,=4 paratodoo ne IN e 
F=*f onde f:A> IR é uma função integrável, então 


“Lo ercgas=[(J, rodar) J(J ro) 


n 


2 


isto é, neste caso, os valores do “integral e do integral usual coincidem. 


Considerando em particular a função interna 4, =[(d, ), E] que, embora te- 
nha sombra nula em !R, não é a função nula em * !R, tem-se o seguinte 


7) = ae poa E = 
Alide= 1 RO) elxJde = (0) 


onde «> é um número real qualquer e q é uma função contínua numa vi- 
zinhança da origem. Estas duas propriedades indicam que a função interna 
A, poderá representar, no universo da ANS, o objecto standard que é a 
função delta de Dirac na medida em que 4, satisfaz algumas das proprie- 
dades mais características de ó, adequadamente interpretadas. Se, no en- 
tanto, se pretender trabalhar com derivadas da função ô então A, (que não 
é “diferenciável) pode não ser a representação mais conveniente. Neste 
caso será desejável representar ô por uma função interna A de tal modo 
que, para todo o me IN,.ô!” possa ser representada por *D“A, Substitui- 
-se deste modo a derivação generalizada de funcionais standard pela deri- 
vação usual de funções internas do universo da ANS. 


3.2. Estudo Elementar de Funções Delta 


Seja g(t), re |k uma função pertencente a C” (IR), positiva” e monotona- 
mente convergente para zero quando ! += tal que 


f eltjdt=1. 


A função interna A :*ik=*Ik definida por A (x)=w.te(wx), onde 


o=[(n), n Je “IN., satisfaz a condição 7 


Aa=0 Woemin (12) 
Dada uma sucessão (u,), de números reais positivos que tende para +=s 
então r=[ta,), 5] é um número hiper-real positivo infinito; como 


1=º elodr=tim [O sle)de=lim P” ng(ur)de 
pode, consequentemente, escrever-se 


if a.c-|(]] neta) 1 (13 


As equações (12) e (13) mostram que a função A, realiza, no universo da 
análise não-standard, a compatibilização das propriedades (1) e (2) referidas 
em $1, após a substituição do sinal = pelo sinal = e do símbolo «s pelo nú- 


* Esta restrição pode ser substituída por outras condições permitindo que g seja ape- 
nas condicionalmente integrável. Deste modo poderão obter-se representações não 
standard para a função delta com sinal arbitrário. 

* A justificação de que esta condição fica assegurada pelas restrições impostas à 
função q, pode fazer-se apelando ao Teorema de Pringsheim mencionado, por 
exemplo, no Curso de Álgebra Superior de J. Vicente Gonçalves (2.º ed. pg 117). 


mero infinito positivo. º Seja v=[(m,)4] um número infinito positivo arbitra- 


riamente fixado (pode supor-se m,>0.Y,.n ); como 


tim [” e(o)dr=1 


então 


1= [Del mar=" [ore oxas = [adaga 


donde se pode concluir que, qualquer que seja o número hiper-real infinito 
v2v/o. 


* evo ke per 
[a aja=* [7 adajár=o. 
Fazendo, por exemplo, v=vw vem v/w=1/Vw=eemon(0), obtendo-se 
então 
se * pe * pe 
dl adode=1 o | adadr=0=" [Pad 


o que mostra que o “integral da função interna A, é infinitesimal fora de 
uma vizinhança infinitesimal da origem. 


Seja agora f:Ik-> uma função limitada º contínua na origem; qualquer 
que seja o número infinitesimal 0<e=[(e, ), SI tal que wee*IR. e qualquer 


que seja o número infinito y> we, vem 
*per Peporporp 
[acrrtga= PPP toe rtagas. 
Como fé limitada então 


ICE ntortga) | < pu l(£ sas) |o 


e, análogamente, para o integral estendido ao intervalo [-y.-£]. Obtém-se 
então 


e per 
A x)* f(x)de 


s * per % 
A da)*f(x)dr= IN Ss dx) f)dr=*f(7) [R A dx)dr=*f(7) 


onde -E<T<+e. Da continuidade de f na origem resulta que se tem 
*f(t)=/(0) e, portanto, 


Se, em particular, se escolher para g uma função indefinidamente diferenciável em 
!k e de suporte contido num intervalo [0, a], para algum a>0, então A (1). xe“ik é 
uma função nula em todo o ponto standard (incluindo a origem) que não é a fun- 
ção nula do universo da análise não-standard. 

No caso de f não ser limitada obter-se-á o mesmo resultado desde que g seja uma 
função de suporte compacto. Neste caso À, será nula (e não apenas infinitesimal) 
fora da mónoda da origem. 
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q 
[st rtadr= 1(0) 


ou seja: A, realiza a operação de amostragem da função f, operação esta, 
que é característica da função delta de Dirac. Consequentemente, pode es- 
crever-se 


<ôp>=a[ | A) *ftajda = 110). 


Analogamente, se obteria a igualdade 


<8",f >=s[ IE *Dºa (x) Ef ve) (EDP fo). (14) 


para me IN, qualquer e f limitada e suficientemente diferenciável na origem. 


Definição 3.1 Chamar-se-á aqui função delta a toda a função interna inde- 
finidamente “diferenciável que satisfaça as propriedades (12) e (13).'º 


Designando por C/[0] o espaço vectorial (sobre « ) de todas as funções li- 
mitadas p-vezes continuamente diferenciáveis numa vizinhança da origem 
então, embora o traço de uma função delta em !R em geral não exista, 
pode falar-se do traço de uma função delta em €,t9) ou do traço da *deri- 
vada de ordem p de uma função delta em c/[0). Os espaços c/[0) fun- 
cionam assim como uma espécie de filtros” do universo standard que per- 
mitem detectar a presença de funções delta ou suas “derivadas no universo 
não-standard onde aquele se encontra mergulhado. Todas as propriedades 
lineares (soma, produto por um numero finito ou *derivação) são detectáveis 
por aqueles “filtros”. O mesmo não se pode dizer quanto às operações não 
lineares. Por exemplo, o produto de funções delta ou de “derivadas de fun- 
ções delta já não é geralmente detectável através dos “iltros” C, (0). 


Exemplo 3.1 Seja 6 uma função positiva indefinidamente diferenciável de 
suporte compacto contido no intervalo [a,b], cujo integral estendido a R é 
igual a 1. Então A,(x)=.*0(mx), xe*IR é uma função delta que se anula 


fora do intervalo infinitesimal [2,.2,]<mon(0), onde e =ajo e e=b/jo. 
Dada uma função /eC,([0) tem-se então 


ii A,(x) Ef(a)dr= (0). 


O quadrado de 4, , é uma função interna indefinidamente *diferenciável que 
em geral já não será uma função delta. Para fec,(0) qualquer 


* por de per HE pur 
) a(o) fed" [Casa *rljde=a,(7) [ado *rtajar= a, (e)r(0) 
*º Note-se que esta definição comporta classes de funções internas mais gerais que 


as apresentadas até agora. Por outro lado, convém notar também que existem na 
literatura definições de função delta ainda mais gerais que esta. 


onde £,<T<e,. Se f(0)*0 então A,(t)f(0) será, geralmente, um nú- 
mero infinito o que mostra que A, não é uma função delta. 


Pode acontecer, no entanto, que O seja uma função para a qual A,(T) 
tome, para cada feC,(0], um valor limitado; neste caso obter-se-á 


%e per * po 
J AG (x) *flx)dr=c,. Ih A,(x) *flx)ds 


igualdade esta que pode servir de justificação para o resultado 5º =cô, com 
c constante arbitrária, que aparece referido na literatura. 


Exemplo 3.2 A função interna definida por 


1 (0) 


A(x) xe*R 


no 


é, como facilmente se verifica, uma função delta. 
O produto xa(x). xe*/k é uma função interna que é infinitesimal em todo o 
*IR, tendo por isso traço nulo em !R ; para toda a função / ec,(0) tem-se que 


(É xA(s) “fas )=0 


o que corresponde ao resultado (standard) xô(x)=0. No entanto o facto de 
xA(x) ser infinitesimal mas não nula é importante em conexão com o Re- 
sultado de Schwartz sobre a impossibilidade de definir um produto de distri- 
buições que prolongue o produto usual de funções. A distribuição vala ai 
valor principal segundo Cauchy de x", pode ser representada pela função 


interna 


mx 


v(xy)= xe *IR 


o” 


tendo-se então que o produto V(x)xA(x) não depende da ordem pela qual é 
executado. Pode, no entanto, provar-se que tal produto, no universo da ana- 
lise não-standard, não representa qualquer distribuição (standard). 


Note-se ainda que 


! 7 se*R 


rox 


difere de 1 por uma função que é infinitesimal em toda a parte excepto na 
mónada da origem; para qualquer peC; tem-se, 


Lo 


e «v(s) “eta =" Ea cota EF *olxjdr *[ 1*o(s)de 
p É 5) ? 


mira 


o que significa que o traço de xv(x) em C; é infinitamente próximo da dis- 
tribuição regular gerada pela função 1. 


Exemplo 3.3 Muitos outros resultados sobre produtos de distribuições que 
aparecem referidos na literatura são susceptíveis de uma interpretação não- 
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standard como a dos exemplos anteriores. Para citar apenas mais um 
exemplo, observe-se que a igualdade 


v(ya(ij=L —et 
T (ros 


o 


está de acordo com o redultado V, [x ']-5 


Note-se porém que não sendo as distribuições representáveis por funções 
internas de forma única as conclusões sobre produtos de funções internas 
não podem transferir-se directamente para o universo standard. Trabalhando 
com representações distintas para V, (x) e à obter-se-ão, em geral, resul- 
tados também diferentes. 
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ESTUDO GEOFÍSICO E 
VIGILÂNCIA DE VULCÕES 
ACTIVOS: O VULCÃO DO FOGO 


Sandra Heleno da Silva 


Aluna do 5.º ano da LEFT e estagiária do Grupo de Fisica da Terra e do Ambiente do IST (GFTA) 


A observação de vulcões em todo o mundo e o estudo das erupções vulcâ- 
nicas levou à descoberta de sinais precursores destes desastres e ao de- 
senvolvimento de métodos de vigilância, que consistem na monitorização de 
um conjunto de parâmetros geofísicos e consequente estabelecimento de 
uma fronteira delimitando a “normalidade” do vulcão. 


Posto isto, a detecção de desvios à “normalidade” (em pelo menos duas va- 
riáveis, para evitar falsos alarmes) deverá levar às acções de protecção civil 
adequadas a cada caso. 


PARÂMETROS GEOFÍSICOS SIGNIFICATIVOS 
Sinais sísmicos 


A sua medição é o método mais comum na previsão de uma erupção vulcâ- 
nica, pois esta é antecedida por um aumento (em número e violência) dos 
tremores locais, devidos quer à ascensão do magma, quer à ocorrência de 
fracturas que lhe dão passagem, podendo mesmo observar-se uma subida 
do foco ou fonte dos sismos — na erupção do Kilauea de Dezembro de 1959, 
os sismógrafos instalados permitiram um aviso com seis meses de antece- 
dência, registando tremores à profundidade de 50 km; com a passagem do 
tempo, a profundidade focal dos sismos foi diminuindo, o que permitiu uma 
medida da velocidade da subida e a previsão da data da erupção. (1) 


Outras evidências precursoras das erupções têm vindo a ser vigiadas e estu- 
dadas, nomeadamente a predominância de ondas sísmicas de alta frequência 
(associada à predominância de formação de novas fracturas em detrimento da 
actividade das já existentes), a diminuição do parâmetro b — declive da curva 
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logaritmica magnitude/freg. de ocorrência (relacionada com o aumento da 
área das fracturas, devido ao aumento do campo de tensões) e o surgimento 
de sinais de frequência confinada e decaimento fraco (que se acredita serem 
gerados por ressonância forçada nas condutas magmáticas). (2, 3) 


Deformação do solo 


Antes de uma erupção é de esperar um entumescimento do aparelho vulcà- 
nico, devido ao armazenamento de magma ou gás; esta deformação pode 
ser medida por inclinómetros, e uma mudança de inclinação que não caia 
dentro das usuais (devidas à variação da temperatura, às marés terres- 
tres,...) deve ser tomada como um sério sinal precursor de erupção, como 
várias observações em vulcões activos têm demonstrado (uma variação gra- 
dual da inclinação do solo precede a erupção, durante a qual se dá um re- 
torno à posição inicial). (4, 5) 

Também a deformação horizontal pode ser monitorizada, mediante a coloca- 
ção de extensómetros em fracturas pré-existentes. 


Outros parâmetros 


A temperatura de nascentes, fumarolas e mesmo do solo sofre muitas vezes um 
aumento brusco antes de uma erupção; este parâmetro deve no entanto ser ana- 
lisado em conjunto com outros, pois dá muitas vezes origem a falsos alarmes. 


Métodos químicos (análise química de fumarolas e nascentes), magnéticos 
e gravimétricos (detectando movimentos do magma em profundidade) têm 
sido usados em menor escala. 


LIMITAÇÕES 


Os sinais precursores de uma erupção variam de vulcão para vulcão : facto- 
res como o grau de obstrução das condutas entre erupções ou teor de sílica 
do magma (intimamente relacionado com a viscosidade das lavas) são de- 
terminantes, tendendo os vulcões de condutas fechadas, alta viscosidade e 
erupções menos frequentes a manifestarem maior sismicidade e deforma- 
ção do solo, pela sua capacidade de suportarem maiores tensões. (6) 


Assim, enquanto num extremo temos, por exemplo, a erupção de Cose- 
guina, que se deu sem qualquer sinal precursor, temos também no outro ex- 
tremo casos como o Vesúvio, que nos últimos anos tem dado fortes sinais 
de aviso sem consequências. (1) 

É portanto de extrema importância um estudo sério dos processos eruptivos 
dos vulcões activos, associado a um esquema de monitorização. Existem, 
no presente, grandes disparidades nesta área para diferentes zonas do 
globo, com uma grande parte dos vulcões potencialmente perigosos sem a 
vigilância adequada. 


O VULCÃO DO FOGO 


Sendo a única a manter um vulcanismo activo entre as ilhas de Cabo 
Verde, a ilha do Fogo é ela própria um enorme vulcão: de forma grosseira- 
mente cónica, descentrada, apresenta no seu topo (a 1600 m de altitude) 
aquilo que é a estrutura mais marcante do Fogo, a Caldeira, uma depressão 
aplanada em hemiciclo, com 9 km de diâmetro, delimitada (excepto a Este) 
por uma escarpa, a Bordeira, que atinge 1000 m de altitude em relação à 
base. Esta estrutura formou-se por colapso do bloco central, e a aproxima- 
ção da superfície à câmara magmática levou ao reinício da actividade vulcà- 
nica, modelando aquilo que na ilha é chamado “o Vulcão”, o cone eruptivo 
principal, que atinge 2829 m. (7,8) 


A base da Caldeira, ou “Chã das Caldeiras”, tem sido o celeiro da ilha , com 
terrenos férteis para a vinha e árvores de fruto, o que levou à fixação das 
pessoas, constituindo as aldeias de Boca Fonte, Portela e Bangaeira. 


Com uma história eruptiva recente (a última erupção tinha ocorrido em 1951) 
e a proximidade das populações, era patente a necessidade de um pro- 
grama de vigilância, e em 1992 deu-se início ao projecto “Estudo geofísico 
do Vulcão do Fogo”, coordenado pelo Prof. João Fonseca, do Dep. Física do 
IST, e consistindo no seguinte: 


* Período de monitorização preliminar (6 meses — 1 ano), com o objectivo 
de definir o patamar de actividade numa situação normal. Após uma insta- 
lação inicial de 3 estações sismográficas cedidas pelo Instituto de Meteo- 
rologia (IM), pelo período de 3 meses, foram montadas 10 estações sis- 
mográficas telemétricas* (cedidas temporariamente pelo Natural 
Environment Research Council do Reino Unido), 4 inclinómetros e 1 ex- 
tensómetro ** (instalados com a colaboração do Institut de Physique du 
Globe de Paris), e ainda 1 sonda geotérmica. ** Esta primeira fase do pro- 
jecto terminou em Agosto de 1994. (9) 


* Monitorização permanente, permitindo a detecção de desvios ao patamar 
definido, e o melhor conhecimento do vulcão. No dia 3 de Abril de 1995 
estavam-se ainda a reunir as condições e os apoios necessários. 


A ERUPÇÃO 


A 3 de Abril de 1995, sem complacências para com os planos humanos, o 
Fogo entra em erupção, pela segunda vez neste século. Tal como aconte- 
cera na primeira, a erupção não se dá na cratera antiga do pico, rompendo 
antes no sopé deste, através de sete focos eruptivos, número que diminui 
com o passar dos dias. A erupção foi antecedida por vários tremores de 
terra, sentidos por alguns dos habitantes da ilha, e apresenta características 
quer explosivas (com projecção de bombas vulcânicas e lapili, localmente 
designado por “jorra"), quer efusivas. 

A erupção não causou vítimas mortais, e apenas algumas pessoas ficaram 
ligeiramente feridas devido ao pânico que se gerou na altura. As lavas, que 
pela sua viscosidade avançam lentamente, alcançaram no dia 8 de Abril a 
aldeia de Boca Fonte, engolindo literalmente as casas, a cooperativa e o 
precioso depósito de água, perante o olhar corajoso dos seus habitantes. 


Destas estações, 6 foram instaladas no Fogo, 3 na Brava e 1 em Santiago. 
Equipamento instalado na Chã das Caldeiras 


19 


20 


A escarpa oeste da Bordeira tinha sido alcançada na noite anterior, cortando 
o acesso a Portela e Bangaeira, para onde se dirigia outra frente. 


No dia em que escrevo, 17 de Maio, Boca Fonte continua a ser a única al- 
deia atingida, e não há razão para pensarmos que o cenário vai mudar. Os 
terrenos agrícolas da Chã foram, na sua grande parte, destruídos. 


O acompanhamento 


A pedido das autoridades caboverdianas, uma equipa do IST (João Fon- 
seca, Sandra Silva e Daniel Lála, Dep. Fisica) desloca-se ao Fogo a 5 de 
Abril, munida de quatro estações sismográficas de “Strong-motion” , cuja 
instalação é concluida a 7 de Abril, nas povoações de Mosteiros, Cova Fi- 
gueira, Patim e S.Filipe, situadas no exterior da caldeira e rodeando esta. 
A escolha dos locais teve de ponderar factores como a geometria da rede, o 
“fácil” acesso (para permitir recolhas diárias dos dados) e a ligação à rede 
eléctrica (que não é frequente na ilha). 

A esta equipa inicial juntou-se a 8 de Abril Lídia Quental, do Instituto Geoló- 
gico Mineiro e estudante de mestrado do Dep. de Minas do IST, e a 11 de 
Abril um elemento do British Geological Survey (BGS), David Stewart, com 
equipamento sismográfico para reforçar a cobertura da ilha — cinco esta- 
ções telemétricas munidas de painéis solares, que foram instaladas em Chã 
das Caldeiras. De 19 a 26 de Abril o projecto contou com a colaboração de 
Sancho Oliveira, aluno do 4.º ano da LEIC e a 27 de Abril chegou ao Fogo 
José Luís Pereira, colaborador do GFTA actualmente a realizar o doutora- 
mento na Universidade de Norwich, com mais 4 estações sismográficas ce- 
didas por aquela instituição. *** 


O objectivo prioritário deste estudo é a protecção civil. As estações sismo- 
gráficas de “strong-motion” permitem a visualização (no campo) dos eventos 
registados, permitindo uma primeira estimativa da localização epicentral (lati- 
tude, longitude) dos sismos, em tempo quase real. Esta informação, combi- 
nada com os histogramas de ocorrência dos sismos é um valioso indicador 
do estado da actividade vulcânica. 


Um bom exemplo da eficácia do esquema montado ocorreu logo no dia 
10 de Abril: a um acréscimo da actividade sismica durante a manhã (culmi- 
nada por dois sismos sentidos pela população em algumas aldeias da zona 
SE da ilha) juntava-se uma forte diminuição da quantidade de lava expelida 
pelo vulcão. Esta situação foi diagnosticada como potencialmente perigosa 
(pois podia significar um forçamento, pelo magma, de uma outra saída), 
pelo que se recomendou ao Presidente da Câmara dos Mosteiros, Hugo Ro- 
drigues, que promovesse medidas adequadas de protecção civil. De ime- 
diato foi ordenado o corte da estrada que leva à Chã e uma hora depois, a 
parede Sul do cone eruptivo ruiu, passando as lavas a ser expelidas nessa 
direcção. A topografia reencaminhou-as para a Chã, e não para a zona ter- 
minal da Bordeira, pelo que felizmente o único resultado deste episódio foi a 
queda de jorra e cinza em toda a zona oriental da ilha. 


No Fogo encontra-se desde o dia 11 de Maio Telmo Silva, Eng. Electrotéc- 
nico recém-licenciado pelo IST, que assegurará a manutenção da rede até 
31 de Julho. A erupção prossegue, com características menos “espectacula- 
res”, pois o magma sai agora por fissuras na base do cone, infiltrando-se 
sob a lava anteriormente expelida; a velocidade da frente de lava estacionou 
nos 2 metros/dia e nenhuma povoação se encontra sob perigo imediato. 


** De 8a 12 de Abril deslocaram-se ao Fogo dois sismologistas da Faculdade de 
Ciências de Lisboa, com equipamento próprio, mas não foi possivel incorporar os 
respectivos dados nesta operação. O mesmo se aplica à presença de dois sis- 
mologistas do IM, que estiveram na ilha entre 7 e 20 de Abril 


Fig. 1 Exemplo de um sismo registado nas estações telemétricas da Chã das Caldeiras — 
Monte Amarelo (FTMA). Monte Cruz (FMCZ), Monte Beco (FMBC), Boca Fonte (FBFT) e 
Monte Losna (FMLN) 


Sugestões para estudos posteriores 


O Vulcão do Fogo apresenta um tipo de vulcanismo ainda pouco estudado 
— o vulcanismo intraplacas, que não pode ser associado nem com as zonas 
de divergência das placas, caso dos Açores, nem com as zonas de conver- 
gência destas, caso do Japão; este tipo de vulcanismo é em geral explicado 
usando a hipótese dos “hot-spots”, pontos em que o manto apresenta “plu- 
mas ascendentes” de material fundido, que “furam” a placa à medida que 
esta se vai deslocando. Mas se esta explicação é intuitiva para arquipélagos 
cujas ilhas são aproximadamente colineares (como o Hawai), já o caso de 
Cabo Verde, em formato de ferradura, não é tão evidente. 


Os dados registados na rede instalada no Fogo vão permitir estudos conducen- 
tes a um melhor conhecimento do vulcão e do seu enquadramento tectónico. 


A distribuição hipocentral dos microssismos registados na Chã das Caldeiras 
durante a erupção dá-nos informação útil no que respeita à localização e 
forma da câmara magmática, partindo do pressuposto que a câmara é uma 
zona de baixa rigidez (portanto assísmica) no interior de uma zona sismica- 
mente activa. (4) 


Se notarmos, entretanto, o facto de a(s) câmara(s) e os canais magmáticos 
serem zonas anómalas de velocidade e atenuação das ondas sismicas, po- 
demos usar uma análise mais sofisticada, denominada tomografia sismica, 
que consiste na modelação da estrutura de velocidade e atenuação por 
baixo do vulcão****. 


A rede de estações vai também permitir o estudo dos mecanismos focais 
dos sismos (orientação do plano de falha, direcção e sentido do desliza- 
mento) e o conhecimento do campo de tensões local (eixos de tensão má- 
ximo e mínimo). 


Podemos assim chegar a um melhor entendimento dos processos eruptivos 
do vulcão do Fogo, o que é extremamente importante para a mitigação do 
risco vulcânico no futuro. 


Assume-se um modelo inicial (parametrizado em células) e calculam-se as pertur- 
bações no modelo a partir das observações do tempo de chegada e amplitude das 
ondas (em várias estações situadas em volta do volume estudado) — como os pa- 
râmetros da fonte sismica não são conhecidos com exactidão, as perturbações nos 
instantes de origem e parâmetros hipocentrais devem também ser calculadas.(10) 
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iti) 


Fig. 2 Resultados preliminares da análise dos dados sismicos registados na rede telemé- 
trica da Chã das Caldeiras; |) localização epicentral dos primeiros oito sismos tectónicos 
analisados: os epicentros estão indicados por pontos, e as estações telemétricas estão re- 
presentadas por asteriscos; indicam-se as localizações da escarpa da bordeira e do Pico, 
bem como o local de erupção (E); ii) projecção dos hipocentros dos oito sismos segundo 
os cortes AA' e BB', com indicação esquemática da topografia; iii) solução de mecanismo 
focal conjunta para OS oito sismos, indicando fortemente a existência de uma falha sub- 
vertical (inclinação de 10") com o azimute de 120" (aproximadamente NW-SE) e actividade 
normal, apontando para a existência de um campo de tensões extensivo na direcção NE- 
-SW e para o afundamento do interior da caldeira. 
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TEORIA DE PERTURBAÇÕES 
E RESSONANCIA - O PENDULO 
SIMPLES 


João Paulo Santos 
Departamento de Matemática do IST 


Nesta pequena nota mostra-se 
como é possivel tratar perturbativa- 
mente o pêndulo simples. Analisa- 
-se a razão do aparecimento de ter- 
mos seculares (termos que crescem 
polinomialmente com o tempo) e 
mostra-se como os evitar. 

Consideremos a equação do pêndulo 


dx 
dr 


==sin(x) 


Para pequenas oscilações podemos 
tentar tratá-la perturbativamente. 
Para tal introduzamos a nova variá- 
vel Y=£+, Desenvolvendo o seno 
em série obtemos a equação 


Ea 
dt 


== Lro(e 
Preto (1) 


Se escrevermos y na forma 


+ex+ 


e introduzirmos na equação (1) ob- 
temos igualando potências de É 


Da primeira equação sai imediata- 
mente que y =Acos(1). Introdu- 
zindo na segunda equação obtemos 


doy 
dr 


4 4 d 
+v =—cos(t)+ — cos(3t) 
x 24 


cuja solução é 


A ã 
w =>>"sin(t)-—= cos(31) 
l6 


192 


Assim o método utilizado não é vá- 
lido pois *, diverge. O problema 


está em que a solução não possui 
periodo 1(a) (ver figura 1). Se es- 
crevermos v(t)=Y(wt) e introduzir- 
mos na equação 1 obtemos 


q yr+Y 


Substituindo agora Y e « na equa- 
ção (1) na forma 


Y 


+EY, + 
wv=|+em,+ 


obtemos o sistema de equações 


Vs tç=0 


Obviamente Y,=Acos(T) (b) com 
t=«w! ficando a segunda equação 


1 


A cos(sr)+| 240, + À lcos(7) 
24 s 


Y+Y 


Escolhendo 


elimina-se o termo ressonante cos(7) 
após o que a equação para Y, fica 


vº+ 


=>-cos(37) 
24 


e obtemos para y 


A 
v=Acos(mi)-£ — 
192 


cos(301) 
com 


4 
m=1- 
16 


Fig. 1 Uma pequena diterença na frequência 
origina um erro significativo para t grande 
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A UTILIZAÇÃO DE POLÍMEROS 
EM FOTÔNICA 


Luís Miguel Lapão 
Grupo de lasers e plasmas 
Centro de Electrodinâmica 


1. AS POTENCIALIDADES DOS 
POLÍMEROS 


Com a diminuição das dimensões 
dos componentes tecnológicos 
cresce o interesse por componentes 
que possam ser feitos à escala mo- 
lecular. 


A óptica não linear é uma nova 
fronteira da ciência e da tecnologia 
que se tem desenvolvido imenso 
com o surgimento das tecnologias 
da fotónica. A fotónica, que utiliza 
os fotões para processamento de 
imagens e de informação, tem sido 
identificada como a tecnologia do 
século XXI, para a qual os proces- 
sos ópticos não lineares proporcio- 
nam as funções chave para a con- 
versão de frequência e para o 
optical switching. 

A fotónica vai substituir complemen- 
tarmente a electrónica. Como são já 
exemplo disso o trabalho desenvol- 
vido em redes de telefones, em dis- 
cos e conecções de computadores 
e em fibras ópticas. 


A este novo conceito correspondem 
vantagens a diversos níveis. Au- 
menta enormemente a capacidade 
tecnológica dos materiais em ter- 
mos de velocidade de propagação e 


de processamento da informação, 
mais, como são fotões que se des- 
locam tem-se portanto a velocidade 
da luz — velocidade máxima de in- 
teracção —, e em termos de densi- 
dade óptica de armazenamento de 
informação. Associados a estas sur- 
preendentes potencialidades co- 
existem a baixa dissipação de ener- 
gia e a baixa interferência crosstalk. 


É no entanto óbvio que a abordagem 
do novo conceito requer uma outra 
forma de pensar e conceber as ar- 
quitecturas dos sistemas, e sobre- 
tudo o repensar dos materiais para 
os novos componentes fotónicos. 


Materiais moleculares e sistemas 
poliméricos em particular, tem sur- 
gido como uma nova e promissora 
classe de materiais ópticos não line- 
ares. Principalmente porque ofere- 
cem propriedades flexíveis, quer ao 
nível molecular quer ainda ao nível 
da estrutura de substância, para a 
optimização da não linearidade e de 
outras propriedades requeridas na 
construção de dispositivos. 


Nos últimos anos, foi possível o de- 
senvolvimento de fontes luminosas 
coerentes de elevada energia, que 
conduziram ao incremento do es- 
tudo sobre as propriedades não li- 
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neares dos polímeros. Os polímeros 
são compostos orgânicos, formados 
por grandes cadeias de monóme- 
ros, fáceis de fabricar e que podem 
ser utilizados, com vantagem, em 
fotónica. 


Estes polímeros devem exibir com- 
portamento funcional relativamente 
à luz: assistidos por campos ou 
mesmo por luz. 


Às grandes vantagens destes mate- 
riais para o desenvolvimento e apli- 
cação do novo conceito advém da 
facilidade em se adaptarem a novos 
arquitecturas, pois são materiais fá- 
ceis de trabalhar, com baixo custo 
efectivo, permitindo uma produção 
em massa. 

Quando se estimulam este elemen- 
tos por meio dos elevados campos, 
fornecidos pelo impulso de um laser, 
as suas propriedades não lineares 
atingem valores mensuráveis a dife- 
rentes níveis. Estes estudos por inte- 
racção laser/polímero pretendem 
fundamentalmente caracterizar as 
susceptibilidades não lineares e pos- 
teriormente, mas com maior ênfase, 
propriedades que permitam o desen- 
volvimento de computadores ópticos. 
As grandes cadeias de polímeros 
dopados que se constroem são ca- 
racterizáveis por: 

* potenciar elevados coeficientes 

não lineares; 


* ter uma grande concentração de 
moléculas activas; 


* permitirem uma rápida resposta a 
campos aplicados (têm portanto 
uma pequena massa electrónica). 


Pode-se pois construir facilmente fil- 
mes finos de multi-camadas, por 
spin coating e por dipping. 

Tem outra propriedade importante 
ao serem compatíveis com muitos 
tipos de substratos: vidro, metais, 
plásticos e silício, HII-IV. 

O facto de terem uma constante di- 
eléctrica substancialmente baixa 
permite-lhes por isso uma resposta 
electrónica muito rápida, em resumo 
apresentam as propriedades de Tai- 
lorhabilidade para a Engenharia mo- 
lecular. Ocorrendo o processo de 
fabrico a temperaturas relativa- 
mente baixas. 


As montagens experimentais con- 
sistem essencialmente em oscilado- 
res ópticos que vão, em condições 
de acoplamento, aumentar a esti- 
mulação das susceptibilidades não 
lineares, tornando assim possível a 
concretização dos efeitos importan- 
tes em fotónica. 


Esta área é obviamente multidisci- 
plinar. Assim o progresso vai de- 
pender da participação activa de fi- 
sicos, químicos, especialistas em 
materiais e de engenheiros. 


RESUMO 


O presente trabalho descreve 
as primeiras fases do pro- 
cesso de selecção e avaliação 
de epífitos criptogâmicos para 
eventual utilização como bio- 
monitores de salinidade at- 
mostérica. Feito o reconheci- 
mento da área de estudo — 
situada no litoral alentejano — 
os líquenes P. periata, R. cali- 
caris é Usnea spp. foram se- 
leccionados para ensaios labo- 
ratoriais: após incubação em 
água do mar artificial, a locali- 
zação celular dos potenciais 
traçadores catiónicos de sal 
marinho (Na e Mgf") foi deter- 
minada por eluição sequencial, 
sendo discutida a correspon- 
dente diferenciação fraccional 
A eficiência na tomada de Na* 
e Cr, em condições reais de 
exposição, veio a ser aferida 
mediante transplantação dos 
líquenes para um local de forte 
influência maritima: a capaci- 
dade de acumulação elemen- 
tal das espécies é discutida 
com base na respectiva morto- 
logia. 


ABSTRACT 


Being a most serious threat to 
materials and/or infrastructu- 
res, sea-salt input over the 
land masses has yet to be gi- 
ven the proper weight it repre- 
sents in the battlefield of air 
pollution: this requires an accu- 
rate, low-cost assessment of 
saline risk at any site, This arti- 
cle presents the first results of 
a research program on biomo- 
nitoring salt deposition density 
by means of cryptogamic 
epiphytes. 
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DE UM PROGRAMA DE 
BIOMONITORIZAÇAÃO 


R. Figueira, C. Branquinho, A.M.G. Pacheco, A.J. Sousa, F. Catarino 


Original submetido em 3 de Março de 1994 


INTRODUÇÃO 


Os líquenes e briófitos são organismos de uma grande simplicidade — 
quando comparados com as plantas vasculares — e têm vindo a ser utiliza- 
dos como monitores de poluição em habitats terrestres (veja-se, por exem- 
plo, [1-10]), uma vez que: 


i) não apresentam cutícula ou camadas cerosas, o que propicia/facilita a 
tomada de água e nutrientes através da respectiva superfície; 

ii) não possuem sistemas radiculares bem diferenciados, o que impede a 
nutrição a partir do substrato, sendo as necessidades metabólicas satis- 
feitas a partir da atmosfera; 

iii) são dotados de forma biológica estável e não-sazonal, de ubiquidade e 
longevidade consideráveis, e mostram ser capazes de acumular elemen- 
tos em excesso dos seus requisitos fisiológicos. 


Neste trabalho, descrevem-se as primeiras fases do estudo que vem sendo 
efectuado sobre epiífitos criptogâmicos, enquanto candidatos a indicadores 
de intrusão salina atmosférica e monitores da respectiva densidade de de- 
posição terrestre: aos materiais e métodos utilizados é feita breve referência 
no número seguinte, passando-se depois à apresentação e discussão do 
percurso inicial de avaliação das espécies consideradas. 


MATERIAL E MÉTODOS 
Área de estudo 


Foi seleccionada como área de estudo e bioamostragem da deposição salina 
uma zona do litoral alentejano relativamente plana e aberta aos ventos ma- 
reiros dominantes, para ocidente da linha Grândola-Santiago-Sines; a exis- 
tência de uma estação de monitorização atmosférica em Santo André pesou 
igualmente na escolha da referida zona. Feito o reconhecimento da área, foi 
assinalada a presença das seguintes espécies: Parmelia caperata (L.) Ach., 
P. perlata (Huds.) Ach., Ramalina calicaris (L.) Fr., R. canariensis Steiner, R. 
farinacea (L.) Ach., R. fastigiata (Pers.) Ach. e Usnea spp. (líquenes); Sclero- 
podium touretii (Brid.) L. Koch (musgo); Frullania dilatata (L.) Dum. (hepática). 


O local-controlo — onde foi recolhido o material para transplante e ensaios 
de laboratório — situou-se na Serra de Grândola, a cerca de 14 Km da linha 
de costa; Ramalina spp. foi colhida em oliveira (Olea europaea), e P. perlata 
e Usnea spp. em sobreiro (Quercus suber). 
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Transplantação de liquenes 


Foi seleccionado como local-piloto para a realização de transplantes o peri- 
metro da estação atmosférica de Monte Velho (Santo André), operada pela 
Direcção Geral do Ambiente e Recursos Naturais (DRARN-Alentejo) e locali- 
zada a cerca de 1100 m da praia. Nos ensaios foram usadas R. calicaris, P. 
perlata e Usnea spp. colhidas no local-controlo, e expostas no local-piloto 
em suporte de madeira a cerca de 1.5 m do solo com orientação de sudo- 
este, Foram recolhidas amostras após um, dois e três meses de exposição. 


O transplante de R. calicaris e Usnea spp. consistiu nos próprios ramos em 
que estes organismos se desenvolviam; P. perlata foi destacada do forófito 
e fixada sobre uma placa de cortiça (20x50 cm) com fio de nylon. Em todo o 
trabalho de colheita, transplantação e processamento laboratorial, foram uti- 
lizados(as) sacos e luvas de polietileno e pinças metálicas. 


Determinação de catiões 


Para a determinação dos catiões, as amostras foram colocadas numa at- 
mosfera saturada em humidade durante 24 h, a fim de serem evitados da- 
nos de membrana por dessecação. Foi realizada a sequência de eluição 
proposta por Brown & Wells [11] para quantificação de fracções catiónicas: 
(1) superficial e intercelular, (2) iões ligados à parede e membrana celulares 
e (3) intracelular. 

Réplicas de 15-20 mg de peso seco (p.s.) foram agitadas duas vezes em 10 
mi de água bidestilada por 40 e 30 min, respectivamente, obtendo-se a frac- 
ção catiónica superficial e intercelular: esta é constituida pelos catiões que 
se encontram depositados à superfície e nos espaços intercelulares do li- 
quene. Em seguida, foram agitadas duas vezes em 10 ml de NiCl, (20 mM) 
durante 40 e 30 min, de modo a eluir por permuta iónica os catiões ligados 
à parede e membrana celulares. Determinou-se o p.s. das réplicas, as quais 
foram depois agitadas em 10 ml de HNO, (1 M) durante 2 h para obter a 
fracção intracelular e, com esta, os iões presentes no interior da célula. 


Os elementos foram doseados por espectrofotometria de absorção atómica 
(emissão de chama) num aparelho Varian Techtron; os resultados totais 
para o sódio resultaram da soma das três fracções e vão expressos em 
umol g' p.s.. 


Determinação de cloreto 


As amostras dos transplantes foram inicialmente colocadas em estufa a 
80 C por 24 h; depois, réplicas de ca. 40 mg p.s. foram agitadas em 10 ml 
de água bidestilada durante 2 h. [CI] foi determinado por titulação mercuri- 
métrica [12], utilizando um titulador automático Metrohm 665 Dosimat; os re- 
sultados vão expressos em umol g' p.s.. 


Incubação salina 


Os líquenes colhidos no local-controlo foram colocados em câmara húmida 
saturada 24 h antes da incubação. Amostras de 100 mg foram agitadas du- 
rante 2 h em erlenmeyers com 1000 ml de diferentes soluções preparadas a 
partir de uma mistura comercial de sais (Sigma), e utilizadas em teores 
equivalentes a 25, 50, 75, 100 e 200% do de uma água do mar típica. Em 
seguida, o excesso de solução presente na superfície das amostras foi reti- 
rado mediante agitação em água bidestilada durante 30 min; réplicas de 15- 
-20 mg p.s. foram então processadas para determinação das fracções catió- 


nicas ligadas à parede celular e intracelulares, segundo o método acima 
descrito. 


Análise estatística 


A existência de valores numericamente anómalos foi pesquisada, com rejei- 
ção sistemática dos outliers de 1.º e 2.º ordem [13]. Foram também conside- 
radas biologicamente espúrias as réplicas com um valor de [K*] no interior 
da célula inferior em 50% ao de controlo, por se considerar que, nessas 
condições, a planta teria sofrido grandes alterações fisiológicas [14]. 


RESULTADOS E DISCUSSÃO 
Reconhecimento de espécies 


A possibilidade de utilizar epifitos criptogâmicos em biomonitorização de- 
pende da sua existência no local a estudar ou da facilidade em efectuar o 
respectivo transplante; as espécies devem ser também de fácil colheita, por 
forma a dispor-se de quantidades suficientes para análise. Com o intuito de 
seleccionar possiveis biomonitores, foi feito o reconhecimento prévio da área 
em estudo, tendo sido identificados alguns macroliquenes e briófitos epifíticos 
(quadro |), depois avaliados de acordo com os seguintes parâmetros: 


A — Presença da espécie na faixa litoral; 

B — Ubiquidade, i. e., distribuição generalizada na àrea de estudo; 
C — Adaptabilidade ao transplante; 

D — Facilidade de colheita e processamento laboratorial; 

E — Capacidade de tomada de catiões (Na”, Mg”, Ca?); 

F — Capacidade de acumulação de aniões (CI). 


O Quadro | mostra que os briófitos S. touretii e F. dilatata têm distribuição 
pouco generalizada na área de estudo, facto que pode ser atribuído à baixa 
pluviosidade na região, com um número de dias biologicamente secos por 
vezes superior a 100 dias/ano [15]; os indivíduos de F. dilatata têm, por outro 
lado, talos muito pequenos e aderentes ao substrato, o que dificulta a recolha 
e limpeza do material: assim, ambas as espécies foram abandonadas no pre- 
sente estudo. A não utilização de P. caperata foi decidida por a espécie reve- 
lar uma fraca capacidade de tomada de catiões, e ser difícil de limpar. 


Quadro | — Avaliação das espécies como biomonitores segundo os parâmetros 
Aa F descritos no texto. (-) ausência; (+) pouco; (++) razoável; (+++) muito. 


Espécies Parâmetros Total 
A B C D E E 
P. caperata - +++ ++ + + ++ 9 
P. perlata + ++ ++ ++ + ++ 12 
Ramalina spp. +++ +++ +++ +++ ++ + 16 
Usnea spp. - +++ +++ +++ +++ ++ 14 
S. touretii - + - + ++ + E 
F. dilatata a + ES ra ++ +» 7 


Foi considerado apenas o género. 


As espécies Ramalina spp., Usnea spp. e P. perlata foram seleccionadas 
para os ensaios laboratoriais e de campo subsequentes, por terem reunido 
o conjunto de características que melhor se adapta à realização de trans- 
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plantes. A. calicaris foi, no género Ramalina, a espécie escolhida por ser 
abundante na área de estudo e de fácil identificação taxonómica. 


Simulação salina 


A selecção de um biomonitor deve basear-se, entre outros parâmetros, na 
capacidade de tomada de um poluente e na resposta diferencial a gradien- 
tes de intensidade e tempo de exposição ao mesmo [5]. Os elementos qui- 
micos podem ter três localizações diferentes nos líquenes [11]. A diferen- 
ciação por fracções — logo, a localização celular — é hoje bem conhecida 
para alguns elementos: o potássio, sob forma solúvel, está quase todo no 
interior da célula e o cálcio restringe-se praticamente à parede celular, en- 
contrando-se o magnésio distribuido por ambos os locais [16]; [K] pode, 
assim, ser um indicador de danos na membrana celular [14]. A determina- 
ção de cada uma das fracções permite, em parte, compreender os meca- 
nismos de sobrevivência dos líquenes face à toxicidade de determinados 
elementos [17]. 


Com o objectivo de avaliar os efeitos da salinidade, os liquenes P. perlata, 
R. calicaris e Usnea spp. foram incubados em diferentes concentrações de 
água do mar. O tempo de incubação (2 h), determinado em experimento 
prévio, foi considerado suficiente para se atingir o equilíbrio entre a concen- 
tração dos elementos em solução e em ligação à parede celular (valores 
não apresentados). A Fig. 1 mostra as fracções extracelulares (elementos li- 
gados ionicamente à parede e membrana celulares) e intracelulares de só- 
dio e magnésio, para cada uma das espécies: pode observar-se que os va- 
lores extracelulares de ambos os elementos são substancialmente 
superiores aos intracelulares. A tomada extracelular parece revelar uma ci- 
nética de saturação para qualquer das espécies, mostrando Usnea spp. 
maior capacidade de permuta para ambos os iões do que R. calicaris e 
P. perlata, com teores de saturação de 100, 75 e 50% em água do mar, 
respectivamente; acima da saturação, Usnea spp. e P. periata apresentam 
uma ligeira diminuição na tomada extracelular de sódio. 


Magnésio 


Na extracelular 


Fig. 1 Teores de sódio e magnésio nos liquenes Usnea spp. (7), R. calicaris (Q) e P. per- 
lata (4), após incubação em água do mar; representa-se a média, o máximo e o minimo 
de três réplicas 


A parede celular dos líquenes possui numerosos lugares de troca iónica com 
carga negativa, aos quais os catiões ficam eficientemente ligados [17]. Con- 
tudo, essa eficiência depende da carga e do raio do catião [18]: Brown & 
Wells [11] indicam mesmo Pb?> Cu?> Ni?”> Zn?> Cd”> Ca”> Mg”> Na”, 
com algumas variações pontuais. A diminuição anteriormente observada nos 
valores extracelulares de Na* em Usnea spp. e P. perlata poderá, assim, fi- 
car a dever-se à competição entre este e outros iões presentes em solução, 
e com maior afinidade para os lugares de permuta. Os níveis constantes de 
K* determinados no interior da célula (valores não apresentados) permitem 
concluir da não ocorrência de danos importantes na membrana celular, o 
que, a acontecer, levaria à redistribuição de elementos pelas fracções celula- 
res e, consequentemente, à perda de significado da eluição sequencial [11]. 


Transplantação de liquenes 


Após a avaliação laboratorial da tomada de traçadores de sal pelas espécies 
P. perlata, R. calicaris e Usnea spp., tornou-se necessário aferir o potencial 
das mesmas em condições reais de exposição: o local escolhido foi, como 
ficou dito, a estação da DRARN-Alentejo em Monte Velho (Santo André), de 
modo a possibilitar futura comparação entre os valores da densidade de de- 
posição salina estimada a partir de biomonitores e por meio de captadores 
atmosféricos clássicos. A amostragem sistemática de líquenes nessa área 
foi condicionada, no entanto, quer por limitação do forófito, quer pela quanti- 
dade disponível dos próprios epífitos; a utilização de transplantes constitui, 
nestes casos, uma alternativa válida, tendo sido já aplicada por alguns auto- 
res na determinação de metais pesados e aniões [19-21]. 
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Fig. 2 Níveis totais de Cl e Na* em transplantes de Usnea spp. (a.b), R. calicaris (c,d) e 
P. periata (e f); representam-se as médias e box-plots das amostras (transplantes: n=30; 
controlo: n=6). 


Os teores totais em Cl e Na* de transplantes colocados em Monte Velho e 
recolhidos ao fim de um, dois e três meses de exposição — durante Junho, 
Julho e Agosto de 1993 — constam da Fig. 2, e são comparáveis aos obtidos 
por vários autores em vegetação natural [22-24]. R. calicaris e Usnea spp. 
revelam boa eficiência na tomada de iões, acumulando diferencialmente ao 
longo do tempo; P. perlata apresenta menor eficiência de acumulação pois, 
ao contrário das primeiras, remete para valores de Cl e Na* que se mantêm 
sensivelmente constantes entre o segundo e o terceiro mês de exposição. 
Este facto poderá dever-se à diferente morfologia evidenciada pelas espé- 
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cies: Takala et al. [24] observaram igualmente valores distintos na acumula- 
ção de CI por Pseudoevernia furfuraceae e Hypogymnia physodes, e refe- 
rem as diferenças morfológicas entre espécies ao comentarem os respecti- 
vos resultados. A menor eficiência de P. perlata poderá ainda decorrer do 
tipo de transplante efectuado, uma vez que esta espécie, dada a estrutura 
foliar do talo, foi fixada sobre uma placa de cortiça e, como tal, exposta uni- 
direccionalmente; pelo contrário, R. calicaris e Usnea spp., que têm estru- 
tura fruticulosa, ficaram suspensas dos ramos originais em exposição omni- 
direccional, o que pode ter favorecido a captação salina. 


Face ao conjunto de resultados obtidos, tanto em laboratório como no campo, 
parece ser inferior a aptidão de P. perlata para a tomada e acumulação de ele- 
mentos traçadores de sal marinho, pelo que esta espécie foi eliminada do 
plano de trabalhos em curso, o qual prossegue com R. calicaris e Usnea spp.. 


CONCLUSÕES 


Os parâmetros que pesaram na eliminação inicial de F. dilatata, P. caperata 
e S. touretii do leque de potenciais biomonitores de salinidade troposférica 
foram, sobretudo, a ubiquidade e a adaptabilidade ao transplante. A simula- 
ção da deposição salina — realizada por incubação em água do mar sintética 
— revelou que a principal forma de tomada de Na* e Mg” pelas espécies P. 
perlata, R. calicaris e Usnea spp. ocorre ao nivel da parede celular, à qual 
os elementos se ligam por um processo de troca iônica. Usnea spp. mostrou 
ser, entre todos os organismos ensaiados, aquele que possui maior capaci- 
dade de permuta para ambos os iões; contudo, mesmo nesta espécie, é 
possível a existência de competição iónica pelos lugares de ligação à pa- 
rede celular. Os transplantes de campo evidenciaram uma apreciável efici- 
ência na acumulação de Cl e Na*, a qual parece ser, no entanto, bastante 
condicionada pela morfologia do líquene. Assim, o inferior desempenho de 
P. perlata, tanto nos ensaios laboratoriais como de campo, leva a que o pro- 
grama de biomonitorização prossiga apenas com R. calicaris e Usnea spp.. 
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RESUMO 


São analisados aspectos im- 
portantes da reacção de oxida- 
ção de o-xileno a anidrido ftá- 
lico sobre catalisadores 
V,OYTIO,, incluindo o meca- 
nismo reaccional e a estrutura 
de sistemas catalíticos. São 
ainda descritos os principais 
processos industriais de produ- 
ção de anidrido ftálico, com 
especial destaque para o 
usado pela GDP - Gás de 
Portugal, S.A. em Cabo Ruivo 
(Lisboa), 


ABSTRACT 


Important features of the oxi- 
dation of o-xylene to phthalic 
anhydride over V,O,TIO, ca- 
talysts, including the mecha- 
nism and the catalyst stucture, 
are analysed. The main indus- 
trial processes for phthalic 
anhydride production are des- 
cribed with particular emphasis 
on the one used by GDP — 
Gás Portugal, S.A. 
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1. INTRODUÇÃO 


O anidrido ftálico é um importante intermediário da Indústria Química apre- 
sentando largas aplicações. Actualmente cerca de 60% da produção mun- 
dial é utilizada na obtenção de ftalatos, que são usados para melhorar as 
características do PVC[1]. É ainda utilizado na preparação de corantes, resi- 
nas, tintas e insecticidas [2,3]. 

O rápido desenvolvimento da indústria de produção de polímeros, verificado nas 
últimas décadas, provocou um forte aumento da procura de anidrido ftálico e con- 
sequentemente da sua produção. Industrialmente, o anidrido ftálico é essencial- 
mente obtido por oxidação de o-xileno sobre catalisadores de pentóxido de vaná- 
dio suportado em dióxido de titânio [1-3]. Presentemente, em Portugal, a produção 
industrial de anidrido ftálico é assegurada pela GDP — Gás Portugal, S.A.. 


2. A OXIDAÇÃO DE O-XILENO A ANIDRIDO FTÁLICO 


2.1. A Reacção Química 


A oxidação de o-xileno a anidrido ftálico é uma reacção complexa que inclui 
a quebra de 6 ligações envolvendo hidrogénio e a incorporação de 3 átomos 
de oxigénio na mólecula de hidrocarboneto reagente (Fig. 1). Além de ani- 
drido ftálico, é referida na literatura a formação de o-tolualdeído, ftalido, dió- 
xido de carbono [1-4] e ainda vestígios de álcool o-metilbenzílico [2], ácido 
o-toluico [2,4], anidrido dimetilmaleico [5,6], ácido benzóico [5,7,8], anidrido 
maleico [6] e monóxido de carbono [9]. No entanto, a nível industrial, são 
obtidos como produtos secundários mais importantes : dióxido e monóxido 
de carbono, anidridos maleico e citracónico e ácido benzóico. 

O complexo mecanismo reaccional de oxidação de o-xileno a anidrido ftálico 
foi estudado por diversos autores [4,5,10-14], apresentando-se na Fig.2 um 
esquema reaccional genérico. Apesar dos numerosos trabalhos publicados 
continua a existir grande controvérsia relativamente à definição das fases da 
reacção e compostos intermediários envolvidos. Assim, alguns investigado- 
res propõem que o anidrido ftálico se forma através de uma série de reac- 
ções consecutivas, envolvendo o-tolualdeido e ftalido como intermediários 
[4,5,14], enquanto que para outros autores [10-13] a sua formação pode 
igualmente ter lugar directamente, num único passo, a partir de o-xileno. 
Para a reacção de oxidação de o-xileno é hoje, de um modo geral, aceite 
que a fase inicial corresponde à adsorção de o-xileno sobre a superfície ca- 
talítica, com abstracção de um dos hidrogénios de um grupo metilo ligado 
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ao anel aromático [1,2]. No entanto, nenhum dos mecanismos publicados 
consegue explicar totalmente a reacção. Nomeadamente, para baixas con- 
versões, é frequentemente referida a formação, sobre o catalisador, de com- 
postos orgânicos pesados, de carácter polimérico, cujas características e 
processos de formação são ainda desconhecidos [4,9,15] 


2.2. Sistemas Catalíticos 


Os primeiros catalisadores usados para a produção de anidrido ftálico a par- 
tir da reacção de oxidação de o-xileno baseavam-se em óxido de vanádio 
puro, os quais eram eficientes apenas para elevados tempos de contacto de 
matéria-prima [2]. Nos últimos anos têm vindo a adquirir maior importância 
catalisadores suportados em material inerte, que lhes confere maior estabi- 
lidade mecância e térmica, aumentando o seu tempo de vida. São referidos 
na literatura diversos suportes dos quais os mais utilizados são a sílica, dió- 
xido de titânio e alumina, podendo a sua estrutura influenciar fortemente o 
comportamento catalítico [16-18]. Actualmente, os catalisadores mais utiliza- 
dos são constituídos por V,O. suportado em dióxido de titânio. Este pode 
existir em três variedades estruturais diferentes, sendo a anatase, fase está- 
vel a baixas temperaturas (T< 973K) [19-22], aquela que proporciona maio- 
res conversões e selectividades em anidrido ftálico [23-25]. A explicação 
para o comportamento deste suporte não está perfeitamente estabelecida. 
Véjux e Courtine [26] explicam o papel benéfico do suporte pelo acordo geo- 
métrico entre planos (010) e (001) da anatase com plano (010) do V;Os ,o 
que favoreceria a orientação dos grupos V=O, tornando-os mais eficientes 
na adsorção e posterior reacção das moléculas reagentes. Esta explicação 
não é, no entanto, unanimemente aceite tendo outros autores associado o 
papel importante da anatase a uma maior dispersão da vanádia sobre a su- 
perfície do suporte [27] e ao efeito estabilizador deste no estado de oxida- 
ção das espécies de vanádio formadas [28]. 


cHa Ea 
+30, — Or& + 3H,0 
cHa “o 


O-XILENO ANIDRIDO FTÁLICO 


Fig.1 Reacção de oxidação de o-xileno a anidrido ftálico 


As espécies de vanádio em catalisadores suportados foram caracterizadas 
recorrendo a diversas técnicas [29-31]. No entanto, apesar dos numerosos 
estudos efectuados não existe ainda um conhecimento total da sua estru- 
tura, a qual parece ser bastante influenciada pelo teor em elemento activo e 
grau de hidratação da superfície. Assim, para teores inferiores a uma mono- 
camada de vanádia superficial verifica-se a existência de espécies isoladas 
de vanádio, as quais podem apresentar estrutura tetraédrica ou octaédrica 
[32-35]. Para teores mais elevados formar-se-iam sucessivamente estruturas 
bidimensionais e estruturas aciculares de propriedades semelhantes ao dió- 
xido de vanádio cristalino [36,37]. De facto, a utilização de espectroscopia 
de raios X [38, 39] mostrou que, para teores elevados em V,Os, não há for- 
mação de camadas sucessivas daquele composto sobre a superfície do su- 
porte. Bond e colaboradores [36] propõem que a formação de V.O. cristalino 
se desenvolveria sob a forma de “torres” ou estruturas aciculares perpendi- 
culares ao plano constituido pela monocamada inicial. 


Os catalisadores industriais são preparados depositando uma fina camada 
do sistema V,O;/TiO,-anatase num material inerte e não poroso (anéis ou 


esferas), o que está directamente ligado ao facto de a reacção ser extrema- 
mente exotérmica, absorvendo o suporte parte do calor desenvolvido [2]. 
A utilização de catalisadores V,O;/TIO, implica, por vezes, a alimentação de 
pequenas quantidades de óxidos de enxofre juntamente com as matérias- 
primas permitindo a obtenção de melhores rendimentos em anidrido ftálico e 
menor desactivação [2]. A presença de SO, poderá favorecer a reoxidação 
do catalisador de acordo com a reacção SO, + V;O, — V.O, + SO,. O teor 
em enxofre a adicionar depende de vários factores nomeadamente do tipo 
de reactor, das condições experimentais e da matéria-prima a utilizar, va- 
riando o teor ideal numa gama bastante larga. 
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Fig. 2 Esquema reaccional para a reacção de oxidação de o-xileno a anidrido ftálico. 


3. PRODUÇÃO INDUSTRIAL DE ANIDRIDO FTÁLICO 


3.1. Processo de Produção 


A produção industrial de anidrido ftálico iniciou-se em 1872, tendo por base a 
oxidação, em fase líquida, de naftaleno com ácido nítrico ou ácido sulfúrico 
na presença de sulfato de mercúrio [3,40]. Este método foi sendo progressi- 
vamente melhorado, passando-se a utilizar, a partir de 1917, a oxidação de 
naftaleno em fase gasosa sobre catalisadores de óxido de vanádio [40]. 


O desenvolvimento da indústria petroquímica proporcionou a obtenção de o- 
xileno em larga escala, a partir da refinação do petróleo, mostrando-se 
aquele composto uma matéria-prima alternativa mais conveniente para a 
produção de anidrido ftálico. Actualmente cerca de 80% da produção mun- 
dial daquele composto é obtida a partir da oxidação de o-xileno em fase ga- 
sosa [1]. De facto, a utilização de o-xileno apresenta diversas vantagens re- 
lativamente à utilização de naftaleno, nomeadamente [1-3,40]: 


— a quantidade de ar necessária para o processo é menor 

— a reacção é menos exotérmica 

— a molécula de o-xileno contém o número exacto de átomos de carbono neces- 
sários não se verificando desperdício de matéria-prima como com o naftaleno 

— pode ser obtido com um grau elevado de pureza a partir da refinação do 
petróleo 

— como o o-xileno se apresenta,à temperatura ambiente, no estado liquido 
é mais fácil de armazenar e transportar. 
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3.2. Produção de Anidrido Ftálico em Portugal 


A produção industrial de anidrido ftálico em Portugal teve início em 1984, 
quando a GDP pôs em serviço a sua unidade fabril localizada em Cabo Ru- 
ivo. Esta instalação, licenciada pela empresa alemã Wacher-Chemie, utiliza 
o-xileno como matéria-prima e tem uma capacidade produtiva anual superior 
a 17500 ton/ano. 


O projecto respectivo, que data do início da década de 80, coincidiu com 
uma importante evolução nos catalisadores à base de titânio-vanádio, que 
passaram a permitir trabalhar com concentrações de o-xileno substancial- 
mente superiores, o que proporcionou grandes melhorias na rentabilidade e 
concepção das instalações. Tal facto foi devidamente tido em consideração, 
pelo que a fábrica portuguesa foi uma das primeiras unidades a usar siste- 
maticamente a relação o-xileno-ar de 60 g/m* PTN em vez dos tradicionais 
44 g/mº PTN. A operação nestas condições, se bem que tenha fortes condi- 
cionalismos de segurança — uma vez que se envereda decididamente pela 
zona de explosividade da mistura — apresenta grandes vantagens económi- 
cas (maior racionalização energética, equipamento mais compacto), a ponto 
de vir a constituir um factor decisivo na sobrevivência das unidades com 
possiblidades técnicas de a adoptar. 


O método de produção de anidrido ftálico utilizado na GDP recorre à oxida- 
ção parcial do o-xileno com o ar atmosférico a cerca de 380ºC e na pre- 
sença de um catalisador muito específico, e que garante uma conversão 
muito próxima dos 78,8%. Esta oxidação, que é realizada num reactor multi- 
tubular (7626 tubos), liberta uma grande quantidade de energia, que é reti- 
rada do sistema através da produção de vapor de água, Uma vez que as 
quantidades em jogo são muito elevadas, a optimização da recuperação e 
posterior utilização desta energia é capital para a economicidade de todo o 
processo produtivo. No caso específico da GDP, a produção de vapor é tal 
que a sua cedência às unidades anexas equilibra praticamente a aquisição 
das restantes utilidades consumidas (energia eléctrica, fueloil, águas). 


Após a oxidação de o-xileno, que é conseguida com uma única passagem 
através do leito catalítico, a mistura gasosa efluente do reactor atravessa 
uma série de arrefecedores onde o anidrido ftálico é separado da multiplici- 
dade de compostos orgânicos co-gerados na reacção. Destes, o principal é 
o anidrido maleico, cuja quantidade produzida (5% do anidrido ftálico) pode, 
em circunstâncias propícias de mercado, justificar a sua recuperação indus- 
trial. Tal não sucede na GDP, onde a totalidade dos sub-produtos orgânicos 
e o CO são destruídos num incinerador catalítico. 

O anidrido ftálico separado em consequência dos processos anteriores é 
posteriormente purificado por dupla destilação, de modo a enquadrar-se nas 
especificações comerciais correntes que estipulam, entre outras proprieda- 
des, uma pureza superior a 99.85 %. 

Ao longo dos 10 anos de funcionamento da instalação verificaram-se algu- 
mas evoluções notáveis nos catalisadores, e de que a GDP tirou o devido 
proveito: desde 1989 que é possível trabalhar com uma concentração de o- 
xileno de 70 g/mº PTN (o que para além de melhorar a rentabilidade do pro- 
cesso, permite aumentar a capacidade produtiva em quase 17%) e, presen- 
temente, está em fase de início de operação um novo tipo que possibilita 
uma taxa de conversão de o-xileno em anidrido ftálico da ordem dos 80.9%. 
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SUMÁRIO 


Populações de partículas exi- 
bindo organização em girofase 
causam acoplamento linear en- 
tre os modos próprios ondulató- 
rios da propagação paralela em 
magnetoplasmas girotópicos. 
Estas interacções permitem (i) 
intensificar instabilidades pree- 
xistentes, (ii) estimular activi- 
dade ondulatória de um tipo 
(por exemplo, electrostática) a 
partir de energia livre (anisotro- 
pia de temperatura) de natu- 
reza distinta (electromagnética), 
e (ii) excitar actividade ondula- 
tória em meios (previamente) 
estáveis. Sendo frequente a 
observação de distribuições de 
partículas não girotrópicas no 
geoplasma, estes efeitos tor- 
nam-se relevantes na interpre- 
tação das sondagens de cam- 
pos de ondas efectuadas por 
missões espaciais. 


ABSTRACT 


Particle populations exhibiting 
organization in gyrophase 
bring about linear wave cou- 
pling among the eigenmodes 
of parallel propagation in gyro- 
tropic magnetopiasmas. These 
interactions can: (i) enhance 
preexisting instabilities, (ij) 
drive wave activity of a given 
type (e.g. electrostatic) with 
free-energy sources of another 
type (e.g. electromagnetic), 
and (ii) stimulate wave activity 
in (otherwise) stable media. 
These effects are relevant in 
the interpretation of space field 
observations because nongyro- 
tropic particle distributions are 
frequentiy encountered in the 
geoplasma. 


SOBRE A ESTABILIDADE 
DE POPULAÇÕES IÔNICAS 
NÃO GIROTRÓPICAS * 


ALL. Brinca ** 
Original submetido em 4 de Maio de 1995 


INTRODUÇÃO 


Populações de partículas não girotrópicas são frequentemente encontradas 
em plasmas espaciais. As suas distribuições de velocidade dependem, no 
plano perpendicular à indução magnética ambiente B,=B,x, da girofase 


o=tan(v./v, ). Electrões e iões com girofases não aleatórias foram detec- 
tados a montante e a jusante do choque em arco terrestre (por exemplo, 
Thomsen et al. [1985], Anderson et al. [1985], Sckopke et al. [1990]), e 
Cairns [1990] observou iões de vapor de água não girotrópicos na vizi- 
nhança do “Vai-vem” espacial. 

Nesta comunicação vamos analisar a estabilidade da propagação paralela 
(k=kx) em magnetoplasmas com uma população iónica não girotrópica, 
centrando o estudo nas frequências abaixo da frequência híbrida inferior. As 
equações de onda e de dispersão mostram que a introdução da não girotro- 
pia pode causar acoplamento entre os três modos próprios de propagação 
paralela em meio girotrópico (ondas electrostática longitudinal, electromagné- 
tica com polarização circular direita e electromagnética ionociclotrónica com 
polarização circular esquerda). Tipos diferentes de não girotropia originam 
classes distintas de acoplamento; em particular, existem organizações de gi- 
rofase que não causam interacção entre os modos próprios girotrópicos, ou 
provocam, apenas, o acoplamento entre os dois modos electromagnéticos 
com polarização circular. Ilustraremos vários efeitos da não girotropia recor- 
rendo a distribuições simples, mas extremas, de girofase (valores iniciais 
idênticos ou diferindo de x). Encontraremos a intensificação de instabilidades 
preexistentes, a estimulação de ondas de um tipo a partir da energia livre as- 
sociada (girotropicamente) a modos de outro tipo, e a desestabilização de 
meios (previamente) passivos. Em todos os casos, a fenomenologia linear 
ocorrida pode ser encarada como resultado do acoplamento entre os modos 
girotrópicos, mediado pelas populações não girotrópicas. 


DISPERSÃO E ACOPLAMENTO 


O magnetoplasma considerado, no seu estado não perturbado, é homogé- 
neo, neutro, e contém uma população de electrões (identificada com o in- 
dice e) e uma população de protões que poderá ser girotrópica (índice p), 


” Comunicação apresentada à Classe de Ciências da Academia de Ciências de Lis- 
boa em 21 Maio e publicada em Memórias da Academia de Ciências, Tomo XXXII, 
ano de 1992/93, págs. 181 a 196. 

” Centro de Electrodinâmica, Instituto Superior Técnico. 
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ou não girotrópica (índice b). As distribuições de velocidades das espécies 
girotrópicas (/ = e, p) são bi-Maxwellianas e definidas por 


Edviv)=M(viA,) 


(mv) (1) 


onde 4=T,,/7,, e Yi =(27,/m) representam, respectivamente, as ani- 
sotropias de temperatura e as velocidades térmicas paralelas (T contém a 
constante de Boltzmann). Magnetoplasmas com A,=1 serão designados 
por isotrópicos, muito embora a propagação não paralela (excluída desta co- 
municação) exiba sempre comportamento anisotrópico. 
Relativamente à população não girotrópica, recordamos [Brinca et al., 1993] 
que a equação de Vlasov para o meio não perturbado homogéneo 
exige distribuições de velocidade da forma F,,=F,(v,.v,.)+82,1), sendo 
2,=q,B,/m, a frequência ciclotrónica da espécie j. Estas distribuições têm 
necessariamente de ser periódicas em f (período 27). Aqui vamos adoptar 


Es = Ev OAB) =20M(v,.A, JD (0 +81) (2a) 


com (& é não negativo e periódico) 


[.oor)dor=1 (2b) 


e definimos os coeficientes de Fourier 


(20) 


Recorrendo às equações linearizadas de Maxwell e Vlasov, uma derivação 
laboriosa registada no Apêndice de Brinca et al [1993] permite, por utiliza- 
ção da equação de onda para propagação paralela, chegar a 


(3) 


onde E é a amplitude complexa (transformada de Fourier em x e t) do 
campo eléctrico ondulatório, os índices + e - identificam as combinações 
0.=[0,ti0.]/2, os elementos da matriz [11,,] valem (/=e.p.b) 


O Dz 
— DD 4 
W-20, kv, (6) (4) 


: (0-29,) 


Si Z(g) 


—0+9,-20, 
ETR 


= E 
> = * Sa 


utilizámos [Fried e Conte, 1961] a função de dispersão do plasma, Z(5), e a 
sua derivada Z'(é), e «, representa a frequência de plasma da espécie j. 

A existência de soluções não triviais para as componentes do campo implica 
a equação de dispersão. 


D(w=,+iy.k)=Det|m, |=0 (5) 


As Equações (3) e (5) evidenciam o acoplamento entre os componentes 
E (0), E(o-2,4),e E (0-22,.k) provocado pela população não gi- 
rotrópica. Esta é responsável pela existência dos elementos não 
diagonais da matriz [m,). A equação de dispersão de um meio girotrópico, 
mmom =0, contém as três equações de dispersão dos modos próprios 
de um magnetoplasma: m,.=0, m,=0,e m =0 representam, respecti- 
vamente, as equações de dispersão dos modos (i) electromagnético com 
polarização circular esquerda para (w.k), (ii) electrostático longitudinal para 


(0-92,,k), e (ii) electromagnético com polarização circular direita para (0-29,). 
TIPOS DE NÃO GIROTROPIA 


A função de distribuição O das girofases da população não girotrópica é 
caracterizada pelos coeficientes &,, do seu desenvolvimento em série com- 
plexa de Fourier (&, é positivo para qualquer distribuição). Os elementos 
da matriz [m,] definidos em (4) mostram que o comportamento linear do 
meio em análise apenas é influenciado pelos coeficientes q. ,e &.,. Caso 
estes coeficientes se anulem, apesar da existência da população não giro- 
trópica (outros coeficientes poderão ser diferentes de zero), o meio exibirá 
propriedades lineares idênticas às do correspondente magnetoplasma giro- 
trópico (girofases aleatórias, uniformemente distribuídas). 

Distribuições de girofase com W,,=0 e 4,40 originam m, =m,=0. Os 
dois modos próprios electromagnéticos não estão directamente acoplados 
entre si, veja-se a Equação (3), mas interactuam através do modo electros- 
tático com o qual estão ambos acoplados. 


Uma situação mais simples ocorre quando os coeficientes de Fourier da pri- 
meira harmónica são nulos (9.,=0), com 4,,40. Verifica-se 
mo =m.=m=m,=0, e o modo electrostático E, fica desacoplado dos 
modos electromagnéticos E, e E (que agora interactuam). 
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Este desacoplamento entre os modos transversais electromagnéticos e o 
modo longitudinal electrostático é fisicamente justificado analisando as den- 
sidades de corrente associadas aos campos transversais. Com efeito, a per- 
turbação linear da densidade de corrente perpendicular (J,., fonte dos 
campos ondulatórios electromagnéticos perpendiculares) apenas pode rece- 
ber uma contribuição da perturbação na densidade (N,, ligada ao campo 
longitudinal electrostático) se a velocidade perpendicular média não pertur- 
bada (<v.. >) for diferente de zero, 


J =eN,<vo >+N <v, > 


Esta condição de interacção implica que os coeficientes de Fourier da pri- 
meira harmónica também sejam diferentes de zero, P. *1), e corresponde 
à existência de uma corrente perpendicular, não perturbada, finita. Plasmas 
girotrópicos, ou populações não girotrópicas com “.,=(), não podem provo- 
car acoplamento linear entre os modos próprios electrostático e electromag- 
néticos da propagação paralela. 

As distribuições de girofase utilizadas nas situações a seguir analisadas cor- 
respondem a 


o alo(o+820)+5(6+82,4+2)] 


0,=Al+e*) (6) 


D=-5(0+02,1). O, (7) 


criando-se, portanto, no primeiro caso, um meio não girotrópico sem acopla- 
mento entre os modos electrostático e electromagnéticos (P, =0). 


INTENSIFICAÇÃO DE INSTABILIDADES 


Começamos por caracterizar o magnetoplasma girotrópico que servirá de re- 
ferência para avaliação dos efeitos da não girotropia. A indução magnética 
estática terá uma amplitude &,=5nT, e as populações electrónica e protó- 


nica, uma densidade comum não perturbada N, =N,=15.5em ': as corres- 


pondentes distribuições Maxwellianas são definidas por 4-1, 
ds 


v,=52x10' kms !,e A Yy =32 kms |. As frequências características 


valem 0, =27x35x10's!, mw, =2mx8I7s' (frequências de plasma), e 
Q =-27x140s', Q2,=27x0.076s'!, e a velocidade de Alfvén toma o valor 
v,=28 kms '. Na apresentação de resultados, normalizaremos as frequên- 


cias (parte real, «,, e parte imaginaria ou taxa de crescimento, y) e núme- 
ros da onda (Kk) relativamente a 82, e Lo, com w,=304kms''. 
v 

O magnetoplasma girotrópico de referência tem energia livre electromagné- 
tica na anisotropia da temperatura dos protões. Esta energia livre estimula o 
crescimento do modo ionociclotrónico com polarização circular esquerda. 
A Figura 1 mostra a dispersão complexa do meio, exibindo soluções numéri- 
cas de baixa (menor do que a híbrida inferior) frequência das equações 


mlok)=m(0k)=m (w,k)=0. De notar que: (i) as duas últimas equa- 


ções de dispersão contêm desvios de frequência de 92, e 282,, respectiva- 
mente, como requerido pela futura análise do comportamento não girotró- 
pico (tem-se 2,=L,, pois a população não girotrópica será de protões); (ii) 


Fig. 1 A dispersão paralela real (1a, painel superior) e imaginária (1b, painel inferior) dos 
modos próprios de baixa frequência, solução das equações m. =m, =m =) parao 
magnetoplasma girotrópico de hidrogénio definido no texto, A energia livre que alimenta a 
instabilidade reside na anisotropia da temperatura dos protões. 


apesar de m (w+20,.k)=D,(w.k)=0 representar, para «, >), a equação 
de dispersão do modo electromagnético (magnetosónico) com polarização 
circular direita, também define, para «o, <0 a dispersão do modo electro- 
magnético (ionociclotrónico) com polarização circular esquerda; e (ili) 
m.(o+82,.k)=0 determina a dispersão complexa dos modos electrostáticos 
ionoacústicos. 

Reorganizando a girofase da população de protões por forma a que a sua 
distribuição passe a ser a indicada em (6), mas mantendo as restantes ca- 
racteristicas do magnetoplasma não perturbado, introduzimos interacção en- 
tre os dois modos electromagnéticos (4...) mas mantemos o modo elec- 
trostático desacoplado daqueles (P, =), A dispersão complexa do meio 
não girotrópico encontra-se representada na Figura 2. Torna-se claro o 


tgp 


Fig. 2 Dispersão complexa do meio não girotrópico obtido do magnetoplasma referenciado 
na Figura 1 através da organização das girofases iniciais dos protões (que assumem ape- 
nas um de dois valores diferindo de 1). 
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efeito do acoplamento entre os modos electromagnéticos: na dispersão real, 
as curvas respectivas acabam por coalescer e assumir um valor comum, 
dm 2, enquanto que a nova evolução das taxas de crescimento conduz a 
uma instabilidade mais intensa e ocorrendo numa banda mais larga de nú- 
meros de onda. 


EXCITAÇÃO DE MODOS ELECTROSTÁTICOS A PARTIR DE ENERGIA 
LIVRE ELECTROMAGNÉTICA 


Partindo do mesmo magnetoplasma girotrópico de referência, com a disper- 
são indicada na Figura 1, mas organizando agora as girofases de maneira a 
que a sua distribuição assuma a forma registada em (7), introduzimos um 
tipo de não girotropia que origina interacção entre os modos electroma- 
gnéticos e electrostático. Na Figura 3 representamos a dispersão complexa 


resultante e a evolução da amplitude das razões |E.(0.4)/E (0-9,.)| e 


um, 


Fig. 3 Caracterização do meio não girotrópico obtido do magnetoplasma referenciado na 
Figura 1 através da organização das girofases iniciais dos protões (que assumem um 
único valor): Dispersão real (painel superior), dispersão imaginária (segundo painel), ra- 
zão [E (0.6)/E (0-2,.4)| (terceiro painel), e razão [E (0-28,.4)/E (0-8,.4)| 
(painel inferior). O modo com frequência real mais elevada (curva a traço-ponto) é instável 
até (valores normalizados) k =40,2, «,=15,5, atingindo a taxa máxima de crescimento 


em 4=30,5, 0, = 


IE (0-20,.8)/E (0-8,.4)|, centrando a atenção apenas nos modos instá- 
veis (com taxas de crescimento positivas). As três ondas instáveis eviden- 
ciam componentes electrostáticos (E) de amplitude relativa importante, em 
contraste com a situação encontrada no meio girotrópico associado, com 
uma instabilidade puramente electromagnética. Constatamos que a introdu- 
ção da não girotropia num magnetoplasma com energia livre electromagné- 
tica (anisotropia de temperaturas) provocou a excitação de modos híbridos 
com intensos componentes electrostáticos. 


DESESTABILIZAÇÃO DE UM MEIO GIROTRÓPICO COM POPULAÇÕES 
ISOTRÓPICAS 


Nos dois exemplos anteriores, o magnetoplasma girotrópico de referência pos- 
suía energia livre electromagnética que estimulava modos ionociclotrónicos. 
Vamos agora tomar como meio de referência um magnetoplasma completa- 
mente estável, sem energia livre, formado por populações (electrónica e protó- 
nica) Maxwellianas e isotrópicas. Trata-se do magnetoplasma definido em 
conjugação com a Figura 1, mas onde se toma 4,=1, eliminando, portanto, a 
anisotropia protónica de temperaturas. A correspondente dispersão complexa 
encontra-se traçada na Figura 4 onde, em contraste com a situação da Figura 
1, todos os modos são estáveis (taxas de crescimento não positivas). 

Neste meio estável, a não girotropia é introduzida criando uma organização 
de girofases com a distribuição (7). Mais uma vez, portanto, se torna possível 
a interacção entre os modos electromagnéticos e electrostático. A Figura 5, 


o 4 8 7 16 20 
Mw Op 


Fig. 4 Como na Figura 1, para o magnetoplasma isotrópico e girotrópico definido no texto. 
O meio é estável 


com o mesmo formato da Figura 3, evidencia os efeitos da não girotropia: o 
meio tornou-se instável, surgindo modos híbridos com taxas de crescimento 
positivas contendo componentes electrostáticos e electromagnéticos. Torna- 


Fig. 5 Como na Figura 3, para o plasma isotrópico e girotrópico da Figura 4, depois da in- 
trodução de organização na girofase inicial dos protões (que assume um único valor). 
O modo com frequência real mais elevada (curva a traço-ponto) atinge um segundo máximo 
da taxa de crescimento em (valores normalizados) k =31. (), =11.1, no valor de 0,362 
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-se, assim, possível estimular actividade ondulatória em meios estáveis atra- 
vés da introdução de não girotropia. 


DISCUSSÃO 


Os exemplos escolhidos ilustram a capacidade das populações não girotró- 
picas para excitar actividade ondulatória inexistente nos correspondentes 
magnetoplasmas girotrópicos. Sublinhamos o carácter linear da fenomenolo- 
gia envolvida nestes processos, o que os torna mais eficientes do que os 
mecanismos de interacção onda-onda susceptíveis de também originar efei- 
tos qualitativamente similares. 

A compreensão do papel desempenhado por populações não girotrópicas 
na estimulação de campos ondulatórios pode justificar a revisão de observa- 
ções efectuadas por missões espaciais no geoplasma. Por exemplo, a de- 
tecção de ondas electrostáticas no pé do choque em arco terrestre [Greens- 
tadt e Mellot, 1987] apresenta dificuldades de interpretação que, 
eventualmente, serão removidas levando em consideração a existência na 
mesma região de populações iónicas não girotrópicas [Burgess, 1987]. 
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RESUMO 


Neste trabalho apresenta-se 
um programa de computação 
para a simulação expedita da 
operação de compressores ge 
neralizado para o caso de mis- 
turas multicomponentes, que 
pode ser utilizado para dimen- 
sionamento preliminar de com- 
pressores ou ainda para revi- 
são de calculo. O programa 
inclui rotinas de calculo de pro- 
priedades criticas das misturas 
gasosas, cujo conhecimento é 
essencial para o dimensiona- 
mento. 


ABSTRACT 


This paper presents a general 
computer program for quick 
and easy simulation of com- 
pressors operation. This pro- 
gram can be used for multi- 
component mixtures, for 
preliminary design of compres- 
sors or tor calculation checks. 
This program includes calcula- 
tion routines for estimation of 
critical properties of mixtures 
which are essential for subse- 
quent design. 
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1. INTRODUÇÃO 


O conhecimento das propriedades criticas de gases é essencial para muitos 
tipos de dimensionamento na Engenharia Química, nomeadamente de uni- 
dades de compressão e de refrigeração. Neste trabalho referem-se alguns 
métodos de estimativa das propriedades criticas reais de gases que foram 
organizadas sob a forma de um programa de calculo automático em lingua- 
gem BASIC por forma a evitar cálculos morosos. Este método é aplicável 
apenas a misturas de hidrocarbonetos gasosos ou a misturas de hidrocarbo- 
netos contendo CO,, H,S, CO e os gases do ar. 

Apesar de já se encontrarem disponíveis diversas correlações e técnicas 
para a estimativa de propriedades criticas de misturas gasosas binarias a 
partir das propriedades dos seus componentes puros, a incerteza aumenta 
quando estes métodos são aplicados a misturas contendo mais do que dois 
componentes. 

As técnicas de estimativa mais vulgarizadas são equações de estado modifi- 
cadas para misturas multicomponentes e o uso do método dos estados cor- 
respondentes, em que se faz a estimativa das constantes pseudo-criticas 
das misturas. 

Contudo, as constantes pseudo-criticas diferem, em geral, consideravelmente 
dos valores reais das propriedades criticas das misturas multicomponentes. 


Note-se porém que, enquanto é possível avaliar os métodos de estimativa 
por comparação com dados experimentais; tal não pode ser feito para as 
condições pseudo-criticas que não tem significado real. 


2. TEMPERATURAS CRITICAS DE MISTURA 
Ao contrário de outras propriedades termodinâmicas, a temperatura critica 
real de mistura não é, geralmente, uma média ponderada com a fracção 


molar da temperatura critica do composto puro, conforme [1]. Li [2] sugeriu 
que se a composição for expressa por : 


9,=),V /z voy 
ata ! (1) 
a temperatura critica real de mistura pode ser dada por : 


1 =Lej, 2) 
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em que: Y, = fracção molar do componente j 
Y. = volume critico do componente j 
T, = temperatura critica do componente j 
T, = temperatura critica real da mistura 
Foram desenvolvidas muitas outras correlações para a estimativa de 7, 


mas esta é uma das mais simples de resolver e é relativamente correcta 
para a maior parte das misturas gasosas [3]. 


3. VOLUMES CRITICOS DE MISTURA 


Muito poucos valores experimentais se encontram determinados, pelo que 
as limitações e a exactidão dos métodos de estimativa ainda não se encon- 
tra correctamente determinada. A correlação de Chueh e Prausnitz [4], mo- 
dificada posteriormente por Schick e Prausnitz [5] parece ser uma das mais 
correctas. Segundo esta, o volume critico é dado por: 


v.=Doy +LLerv, (3) 


com: 

9,=y) (Em, (4) 

V, é um parâmetro de interacção tal que v,=0 e V,, com ;zj é dado 
pelas seguintes equações: 

+Bô, +06) +D6'+E6! (5) 

Pu ffv,+v) (8) 

8,=(v,*-v, AA =V,*4) (7) 


Os cinco coeficientes referidos na primeira destas equações encontram-se 
tabelados [6] para algumas misturas binarias e apenas quando O < à, < 0,5. 
Mesmo para muitas misturas binarias estes coeficientes são iguais a zero e 
não foram ainda propostos coeficientes para misturas com mais do que dois 
componentes. Neste trabalho consideraremos que V, com (i+;)=0. 


Foram testados alguns casos para misturas binarias tendo-se concluído que 
esta técnica introduz um erro sempre inferior a 5%. 


4. PRESSÕES CRITICAS DE MISTURA 


A dependência das propriedades criticas de misturas da fracção molar é 
muitas vezes não linear e a estimativa é, por vezes, de pouca confiança, 
como fizeram notar Reid et al. [1]. 

Uma das melhores correlações parece ter sido a desenvolvida por Chueh e 
Prausnitz [4] que é baseada numa equação de estado tipo Redlich-Kwong 
modificada: 


po=(RT/v.-b)-(a,/riv(v. +b,)) (8) 


sendo T, e V, calculados pelos métodos descritos anteriormente. 


Os coeficientes de mistura para a determinação de ?, são definidos como 


segue: 
b=59b=D(y,b RT./P,) (9) 
a,= 5 E p3, 4 (10) 

com: 

Q, =0,0867-0,01250,+0,011 0) (11) 
a= DRT, (12) 
a,=(2, +82, )RT, (V, +V. )/2(0,291-0,04(09,+00,)) (13) 
To=(1-k MT, T.) (14) 
o, =(RT, Jv.-b)-r)pv.(v +b)/(em, ) (15) 


O parâmetro de interacção 4, varia de 0,1 a 0,01 mas só se encontra tabe- 
lado para misturas binarias. Aqui tomamos um valor médio de &,=0.055, 
Spencer et al. [3] chamaram a atenção para a existência de um desvio mé- 
dio entre ?. e valores experimentais de 2 atm o que é um valor muito redu- 
zido para erro, 


5. PROGRAMA DE CALCULO DE PROPRIEDADES CRITICAS DE MIS- 
TURAS 


Pelas equações atrás apresentadas torna-se evidente que a estimativa de 
propriedades criticas de misturas multicomponentes envolve cálculos muito 
morosos, especialmente quando do calculo dos parâmetros de interacção. 
Assim, organizaram-se estas equações sob a forma de um programa de 
computação que foi testado para duas misturas gasosas binarias e os resul- 
tados obtidos comparados com valores experimentais como segue : 


a) mistura binaria de etano e benzeno com 39,2% de etano. Para esta mis- 
tura tem-se [7]: P =83,8am e T =225,9ºC. A partir das propriedades 
criticas dos componentes puros [1] e a utilização do programa de calculo 
atrás referido, obteve-se: 


P,=76.2atm, T, =219,8ºC é V, =224,9 g/cm' 


Os erros de estimativa são de -2,7% para a temperatura e de -9,1% para a 
pressão. 


b) mistura binaria standard de ar (21% O, e 79% N,) para o que se tem: 
P,=37,2atm e T,=132,3K . Pela aplicação do programa vem: P, =42,7 atm; 


31,3K e V, =86,5g/cm', a que correspondem erros de estimativa 
de 14% para a pressão e de -0,75% para a temperatura. 
O programa irá ainda ser testado para misturas gasosas com mais do que 


dois componentes, quando existam valores experimentais das propriedades 
criticas das misturas para comparação. 
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Entretanto, parece tratar-se de um método útil para a estimativa rápida de 
propriedades criticas de misturas gasosas e julga-se que o erro na estima- 
tiva da pressão critica para misturas binarias não deverá exceder 15%, es- 
pecialmente quando se considerem misturas de hidrocarbonetos. 


6. PROGRAMA DE CALCULO DOS PARÂMETROS DE FUNCIONA- 
MENTO DO COMPRESSOR : 


O programa desenvolvido baseia-se no algoritmo de calculo de Bogart [9] e 
nos métodos de análise de eficiência de Urick e Odom [10] 

A eficiência do compressor é dada por uma série de parâmetros como o 
caudal volumico: 


Lj) 
M 


(16) 


e o volume especifico V pode ser obtido a partir de correlações generaliza- 
das do factor de compressibilidade Z: 


V=R'TZIP (17) 


sendo R' a razão entre a constante dos gases perfeitos e o peso molecular 
médio. 

O factor de compressibilidade é calculado como habitualmente a partir de 
dados tabelados, em função da pressão e da temperatura critica, ou ainda, 
conforme Urick e Odom [10], por correlações envolvendo a temperatura e a 
densidade reduzidas, por um processo iterativo: 


Z=1+6 (bh +h,/T,+b/T)+ Cb, +b/T, +b/T;) (18) 


G=PJT, 2 (19) 


O grau de convergência é fixado pela introdução de um parâmetro -Z. 

Por outro lado, a pressão criada no compressor H é uma função de 

a) razão de compressão: r=Pd/Ps (20) 

b) coeficiente politropico: v=(k-1)/k (21) 
com: k=C,/€ 

c) temperatura absoluta: T 

d) eficiência de compressão: h; 


na seguinte forma: 


H=2RT(o-W/(vn) (22) 


em que 
dia (23) 


Podem utilizar-se vários tipos de eficiências para avaliar a eficiência de 
compressão, tais como a eficiência politropica e a adiabatica. relacionadas 
entre si por 


Na =(r"=)/(r" 1) (24) 


O trabalho necessário para a compressão do gás é dado por: 


W=FH/n (25) 


Deve notar-se que muitas das variáveis envolvidas são dependentes da tem- 
peratura. Embora se possa efectuar o calculo para acréscimos diferenciais 
de temperatura, convencionalmente a equação (22) é resolvida para a tem- 
peratura de sucção Ts e a média aritmética dos factores de compressibili- 
dade para a descarga. Odon [10] recomenda que seja calculado a tempera- 
tura média entre a descarga e a sucção. 


A temperatura de descarga é dada por: 
Td=q Ts (26) 
Calcula-se ainda o número de andares de compressão por : 


1=HIh, (27) 


sendo 4, a perda de carga a vencer por andar, valor standard dado pelo fa- 
bricante. 


Este programa integra assim o programa anteriormente referido [8] para a 
estimativa de propriedades criticas e o algoritmo de calculo de Bogart [9]. 
Tanto quanto temos conhecimento este último algoritmo foi utilizado no cal- 
culo de compressores de amoníaco puro mas não para misturas gasosas 
multicomponentes. Para esse caso particular, e para a compressão de 
100000 kg/h de amoniaco de 265 a 1300 KPa obtiveram-se diferenças de 
estimativa dos parâmetros do processo nunca superiores a 1,5% quando 
comparados com os resultados obtidos pela aplicação do método convencio- 
nal de Mollier [11] à mesma situação. 


O programa global foi utilizado no dimensionamento preliminar de 2 corn- 
pressores, um de gases residuais do fabrico de estireno a partir de etilben- 
zeno, contendo dioxido de carbono hidrogénio e hidrocarbonetos de 1 até 
3,5 bar. 


Foi ainda utilizado no calculo de um compressor de gás combustível com- 
bustível obtido no processo de carbonização da madeira [12], à base de mo- 
noxido e dioxido de carbono, azoto e hidrocarbonetos de baixo peso mole- 
cular, de 1 até 10 atm. 


Futuramente o programa irá ser testado na simulação da operação destes 
compressores. e poder-se-à assim avaliar quantitativamente da sua preci- 
são. 
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ABSTRACT 


In this paper a brief discussion 
of the theory of the simple ro 
tating disc electrode is presen- 
ted. The great advantage of 
the rotating disc electrode Is 
that in many cases one can 
calculate exact solutions to 
many problems in electroche- 
mistty, reaction kinetics, and 
analysis. 
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1. INTRODUCTION 


In the study of the mechanism and kinetics of electrode reactions it is neces- 
sary to make use of a system where the transport of species to and from the 
electrode may be controlled and varied in a known way. Methods involving 
stationary electrodes in contact with an unstirred solution generally give irre- 
producible steady state results and seem to be governed by Albery's [1] 
equation. Thus, to obtain reproducible results, it is essential to employ a 
system where the mass transport of electroactive species to the electrode is 
well defined. 


Several electrode systems have been devised to achieve controlled and pre- 
dictable mass transfer. There are systems where the electrode is in motion, 
for example: rotating disc, rotating wires, dropping mercury, etc.; and others 
where the solution is forced past the stationary electrode: wall jet, bubbling 
electrodes, etc. In this paper, the rotating electrode geometry was preferred 
in the form of a simple rotating disc electrode (RDE) as shown in Figure 1. 
The supply of electroactive species to this electrode is amenable to rigorous 
theoretical treatment and a brief discussion of this theory is presented below. 


J Insulating 
mantle 


A Brass 


Elelectrode material 


Fig. 1 The rotating disc electrode 


2. THE GENERAL MASS TRANSFER EQUATION 


The general rate of mass transfer of species A to an electrode in solution is 
given by the Nernst-Planck equation: 


J=ev-Degrade (1) 


where 7 is the flux of species A (mol S- cm?), D is the diffusion coefficient 
(cm? s7), c and z are the concentration (mol cm?) and charge of A respecti- 
vely, q is the potential and v is the velocity (cm s!). F, R and T have the 
usual meaning. The three terms on the right-hand side represent the compo- 
nents of the flux due to convection, diffusion, and migration respectively. For 
solutions containing an excess of background electrolyte the migration term 


can be neglected. 


3. VELOCITY PROFILE AT AN ARDE 


In order to be able to solve the general equation (1) explicitly, one must first 
know the functional dependence of the fluid velocity v near an RDE. This 
problem was solved approximately by Von Karman [2], and the solution was 
improved by Cochran [3]. Because the rotating disc system possesses cylin- 
drical symmetry, it is convenient to choose cylindrical polar coordinates, r, q 
and x as shown in Figure 2. The corresponding fluid velocity components, 


symbolized v,, v, and v, are given by: 


rw F(X) (2) 

Vo SEWG(N) (3) 

vo==(w uv) HIX) (4) 
(o) 


Vx 


jf 


Fig. 2 Cylindrical polar coordinates and fluid velocity components for the RDE 


where v is in cm s”, w is the rotating speed in radians s!, v = 2 is the ki- 
p 


nematic viscosity in cm? s!, X = «and X,= (2) : X, is the hydrodynamic 


boundary layer thickness (typically 0.1-1mm) and is characteristic of the system. 
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Notice that v, is negative and independent of r. This means that the flow is 
towards the disc rather than away from it, and that the surface of the disc is 
uniformly accessible to material transported from the solution. 


Using the functions F(X), G(X) and H(X), depicted in Figure 3, we can des- 
cribe the fluid flow below the rotating disc. G(X) describes the spinning mo- 


10 


0.5 


0 10 20 ao “0 
X =x/Xy 


Fig. 3 The functions F(X), G(X) and H(X) which describe the fluid flow at the RDE 


tion and falls steadily as X increases. At the electrode surface, G(X) is unity 
and F(X) and H(X) zero. This means that there is a stagnant layer which ro- 
tates with the disc. The radial flow, F(X), has a maximum close to X=1 and 
then decreases further away from the electrode. The spinning disc drags the 
fluid near its surface along with it and, because of centrifugal force, flings the 
solution outwards from the centre in a radial direction. As pointed out previ- 
ously, one very important feature is that v, is independent of r, so in any 
plane parallel to the disc surface v, is constant. It is this feature of the 
hydrodynamics of the RDE system that makes the mathematical description 
of the transport to the disc much easier than that of other forced flow 
systems. The different aspects of the flow pattern, qualitatively described 
above, are depicted in Figure 4. 


4. THE MATHEMATICAL SOLUTION 


From equation (1), after the simplification due to the addition of excess back- 
ground electrolyte, and taking into consideration the continuity equation, we 
can write 


) 
Ce pdivgrade-T. grad « (5) 


dr 


This equation written in cylindrical polar coordinates is given by 


de dp) de 19 * de de de 
=D( +— + 5—5)-v -Vy =v (6) 
dr dr EE E dp” dr dq dx 


Levich [4] used the velocity profile determined by Cochran and solved the 
mass transport equation (6) for the first time. 
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Fig. 4 Schematic representation of various aspects of the fluid flow pattern for an RDE 


Equation (6) is greatly simplified after considering the following boundary and 
symmetry conditions: 
c 


(i) Under steady state conditions de = 
dr 


(ii) For reasons of symmetry, derivatives with respect to º may be ignored. 


(ii) É o atr= 0. Otherwise the Lã term goes to infinity as 1>0. 
Er ER 


7) 

As v, is independent of r, e remains zero until the edge of the disc is rea- 
ER 

ched. Equation (6) is then reduced to 


is [) 


With the condition that far away from the electrode c >c., (where x-00), 
it can be shown that 


(8) 


where X, is the diffusion layer thickness in cm: 


X,=0.643 Wo vit pr (9) 


and W is the rotation speed in Hz. 

Equation (8) relates the concentration gradient at the surface of the elec- 
trode to the surface concentration c,, the bulk concentration c. and the 
length Xp. Figure 5 shows a typical plot of c as a function of distance. It can 
be seen how the length X, divides the concentration profile into two regions. 
For x < X, transport by diffusion is dominant. For x > Xp the solution is well 
stirred and the concentration is uniform and equal to the bulk value c.. 
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Fig. 5 The variation of concentration with distance from an RDE 


5. THE LEVICH EQUATION 
When a species is reduced on the electrode in the steady state, the flux of 
material reacting at the electrode surface must equal the flux being transpor- 
ted through the diffusion layer and arriving on the electrode. For the elec- 
trode reaction we have 

j=k o, (10) 
where k', the rate constant, is given by [5] 


k = Keo exp (-oFE/RT) (11) 


Ke-o is the standard electrochemical rate constant (cm s'') and « is the elec- 
tron transfer coefficient. For the transport, from equation (8) we have 


=€) (12) 


(13) 


where k, is the diffusional rate constant given by 


kp=5— (4) 


Equation (13) has two limiting forms: 


() 1 k5>>K,j=K c.. The flux is under kinetic control. 


(ii) If ko <<k,j=k c.. The flux is controlled by mass transport. (The kine- 
tics are so fast that every oxidized species reaching the electrode is im- 
mediately reduced, so c, = 0). 


Figure 6 shows the current voltage curve for the reduction of an electroactive 
species at an RDE. The two limiting cases are clearly seen. 


Mass 
transport 
control 


—E—— 


Fig. 6 Schematic polarogram for the reduction of an electroactive species at an RDE. Con- 
centration profiles are also shown 


Under conditions of mass transport control (k' >> ks) the limiting flux j, is gi- 
ven by 


ji=kpe =D (15) 


using the fact that j = i/nFA, where n is the number of electrons transferred 
and A is the electrode area (cm?), and equation (9), we can write an expres- 
sion for the diffusion limited current, i,, which is the well-known Levich equa- 
tion [6]: 


à, =1,554 nFAD” ve W" (16) 


This equation provides an excellent test of whether or not the current is enti- 
rely mass transport controlled. If it is, the current will be proportional to W'?. 
In conclusion, it has been shown that the rotating disc electrode has control- 
led and calculable mass transport. This means that we are able to calculate 
the current for the electrode under Kkinetic control, the transport limited case 
and all the intermediate cases. 


6. NON-STEADY-STATE TECHNIQUES 


Although the main application of the RDE is to study the kinetics of a reaction 
under steady-state conditions, various other techniques have been used where 
the concentration profile and current are functions of time. Some techniques in- 
volve the observation of the disc current following a step on the disc potential, 
others consist of the modulation of the disc current, disc potential or even rota- 
tion speed and the comparison of the response to the applied stimulus. 

The solution of the non-steady-state problem for a potential step at the ADE 
has been accomplished by approximation methods [7, 8] and by digital simu- 
lation [9]. 

Finally, hydrodynamically modulated techniques have been described in the 
literature [10-12]. The sinusoidal hydrodynamically modulated technique, for 
example, consists of a small modulation about a central rotation speed. This 
technique is used in a number of ways such as: the determination of very 


59 


60 


low concentrations, studies in the presence of surface complications, and for 
measurements near the background limits for oxidation or reduction. 

Some of these techniques have been used in several studies and more de- 
tails, namely on their practical applications, will be given in subsequent papers. 
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